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ABSTRACT 


A portable three-dimensional computer graphics 
software package was developed utilizing the Fortran 
language. The package included the capability of 
displaying any object as a wire-frame image, as a 
wire-frame image with hidden lines removed, or as a 
solid figure with hidden surfaces removed. This 
computer graphics package provides the user with the 
ability to rotate, scale, and translate any part of 
the displayable image. It was utilized to display 
images or four distinct display devices with only 
minor software alterations. Three totally different 
host computers supported the four display devices. 
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X. INTRODU CTI ON 


Presentation of information using computer aided video 
graphical oi pen plotting devices has become very important 
and extremely useful in almost every profession. 
Three-dimensional computer graphics is being utilized to 
display air combat training in the military, xray scans of 
the human body in medecine, and blue prints and stress 
characteristics for mechanical parts in industry. It has 
even been used to make animated movies. Since one 
manufacturer of graphical machines has net cornered the 
market, software and hardware standards have not been 
established. Each machine has a different screen area, 
smallest resolution size, and data structure. The lack of a 
standard graphics language becomes a source of expensive 
software re-writing every time a newer more capable graphics 
machine is purchased. 

The existing video graphical devices can be divided into 
these four distinct categories: direct view storage tubes; 
vector generator cathode ray tubes (CRT's); raster scan 
CRT's; and plasma panels. The first three types of devices 
have acheived high resolution displays and are most commonly 
used. As a minimum, the hard copy devices have the 
capability of pen plotting and the more advanced machines 
produce shaded images using electrostatic plotting. Except 
fer a plasma panel device, access to the above graphic 
machines was readily available at the Naval Postgraduate 
School. Although all of these devices were supported 
through Ecrtran IV, the software for one machine exhibited 
little discernible similarity to that for another. Each 
graphics device did have its own manual, but the actual 
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efficient utilization of an output terminal, at or near its 
designed capability, was typically left for experimental 
realization. 

The experimental process, learning one machine's 
software idiosyncrasies, required an amount of time which 
was generally not reduced when learning these of a second. 
The high cost of software and this lengthy device learning 
process were the two main incentives urging the development 
cf a graphics software package which could he used with any 
display device. Since Fortran IV has become a universally 
supported language, it was selected as the development 
software. Graphical presentations normally involve one of 
the following four types: two-dimensicnal (2-D) graphs; 2-D 
images; three-dimensional (3-D) graphs; or 3-D images. The 
first two types of presentations have been well documented 
and the theory for both has been fully developed. The 
theory ter and usage of 3-D graphs has also received 
considerable attention. Programs for their display have 
been published in many software languages. Additionally, 
the Naval Station in Keyport, Washington was interested in a 
real time 3-D Torpedo presentation of the torpedo test area 
for the Bange Safety Officer. 

For these reasons, the scope of this research was 
limited to developing a three-dimensional (image) graphics 
package, written in Fortran, which would develop a data set 
that could be displayed on any selected device. A portable 
graphics package would reduce software costs for interfacing 
with a new device to a minimal effort. Usage of a 3-D 
graphics software throughout a large organization, like the 
Navy, would reduce a programmer's learning experience to a 
cne time effort. An individual would then become a portable 



II. THEEE-D IHENSION AL GRAPHICS 


The proper display of a 3-D object on a 2-D surface, 
such as a video screen or a piece of paper, ty a computer 
required that a complete numerical description of the 
object's boundaries or surfaces be supplied for processing. 
The two generally accepted methods used to represent 3-D 
objects are: 

1. " surface definition using mathematical eguations; 

2. and surface approximation by planar polygonal 
mosaic." ^ 

Either description required that a coordinate system be 
constructed to provide the numerical values. 


A. DISPLAY CF AN OBJECT AND ITS HOTION 


A right-handed, cartesian coordinate system provided an 
acceptable and most generally understood system to describe 
3-D objects. The unit of measurement was aribitrary with 
the only requirement that it was consistent. This system 
has been labeled the Object coordinate system. The 
complexity of the irathematics required to define the 
surfaces of an object with equations was only surpassed by 
the algorithms required to project the image onto a 2-D 
surface or alter the viewing aspect. It was much easier, 
both for understanding and utilization, to consider each 
object as a set of one or more polyhedra. Since a 
polyhedron consists cf four or more planar sides, called 
polygons, the realistic approximation of non-planar 
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surfaces, such as a sphere, was determined ty the number of 
polygons used. By definition, a polygon was a flat surface, 
a plane, hounded by three or more connected, straight lines, 
called edges. Two edges intersected at a single point, a 
vertex of the polygon. 

Using the object coordinate system, a vertex was defined 
by its x, y, and z values, with the origin placed at the 
geometric center of the object being described, the location 
of the vertices was simplified. By indexing each vertex as 
it was specified, an edge was determined by the indices of 
the two vertices which were its endpoints. Similarly, 
indexing each edge allowed the description of a polygon 
using the indices of the edges forming its boundary, 
finally, a polyhedron was specified with the indices of the 
polygons which formed its planar surfaces. 

With this type of object definition, object motior was 
inplemented by moving its vertices in a linear manner, 
linear movement of the vertices preserved the straightness 
of the polygonal edges and the structural purity of each 
polyhedron. Translation, which was the linear displacement 
of an object, was defined by the subtraction of a distance, 

I , T , or 1 , from the respective vertex coordinate value, 
x y z 

x, y, or z. Thus, translation produced an entirely new set 
of vertex values, which could be stated as: 



z* = z - T 

z 


Translation coupled with rotation of a point provided 
the ability to describe any 3-D motion. Rotation of a point 
P was most easily explained in two dimensions. Thus, in 
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Figure 1, F was rotated through the angle theta about the 
origin into the point P' by the transformation: 

x' = x cos © * y sin © 

y' = x sin 0 - y cos © 

This derivation for 2-D rotation, as provided in Figure 1, 
was directly applied to 3-D rotation about the Z-axis. The 
angle of rotation was measured in the clockwise direction 
looking from the positive infinity of the axis about which a 
point was to be rotated. Confining the rotation of points, 
and thus objects, to a combination of 2-D rotations greatly 
simplified the computer implimentation. 

The last transformation utilized to alter the viewing 
aspect of an object’s image was scaling. This algorithm 

required that the vertex coordinate values be multiplied by 

the scale factors 3 , S , and S . Frcvided the scale 

x y z 

factors were of equal magnitude, the scaling was linear and 
preserved the polyhedron's shape. 

To project a 3-D object onto a 2-D surface required two 
transformations. The first, called the viewing 

transformation, mapped the object coordinates into a system 
which had its origin at the viewpoint* or the eye, cf the 
graphic software's user. This coordinate system preserved 
the object's linearity and produced the image of the object 
as seen by the "eye". Hence, it was called the Eye 
coordinate system. Its z -axis was used to represent, or 

measure, the depth of the images. The system's X and y 

e e 

axes were aligned with the horizontal and vertical 

dimensions of the display screen, respectively. As shewn in 
Eigure 2, the viewing transformation constructed a 
left-handed cartesian coordinate system. 
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The second transformation simply projected the eye 
coordinate points onto the plane of the display screen. 
This entire transformation was easily constructed and 
explained in Section III. E. 2., using Figure 12. 
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6+<fi+P =180° 

COS/? = COS Q) S\Up- Sw(7T-0-0) 

y'M = 'COS (<£ +Q) % = Sim C0 + 6) 

= -COS0COS0 + SIK/0SW# = COS(/)S\N0 + Sl*«/>COS0 

Y /6 = ~%? cos# + % = % Siiy/0 t % COS0 

v /= xsiin/ 0 -ycos<9 X'= y cos^ + Vs/ a/0 


Figure 1 - TWO-D ROTATION OF A FOINT 
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Figure 2 - EYE COORDINATE SYSTEM 
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IMAGE REALISM 


The perspective projection of an object onto a screen 
often produced points, and thus vectors, which could net be 
displayed. These vectors could be located behind the 
observer cr simply off the screen. The first case would 
have caused erroneous vectors to be displayed, which cculd 
not actually te seen. Typically, the second would generate 
program failure, because the machine could not display 
points, cr vectors located off the screen. Therefore, the 
non-viewatle portion of any image had to te eliminated or 
cut away from the viewable section. This procedure was 
called image clipping. 

Inital 3-D displays portrayed objects as wire-framed 
images. A polyhedron presented in Figure 3, was shown with 
all edges displayed. Because this was a simple body, the 
dedicated observer usually recognized it as a representation 
cf a 3-D oDject. However, the correct viewirg aspect (i.e. 
which surfaces were closest to the viewer) could not 
positively be ascertained. For this reason, the first 
effort to have a computer determine a more realistic 
presentation was the hidden line removal algorithm. Using 
one of these algorithms, the display shown in Figure 4 was 
drawn. The proper viewing aspect was instantly apparent 
provided the viewer recognized the display as a 3-D image. 

The next computer graphics effort to increase display 
realism was to have the computer generate solid polygonal 
surfaces. This required that the display device had a 
shading or full color capability. For realistic images, the 
computer had to possess either a software or hardware 
implemented hidden surface removal algorithm. Now, the 
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viewable image's surfaces were displayed as shaded or 
colored polygonal planes with user supplied shading or 
cclcrs. 

Host recent reasearch has concentrated on producing 
realistic, computer generated shading algorithms. Using the 
most complex shading procedures particular elements of image 
realism have all been acheived. However, a single universal 
solution to the proper shading of images has not been 
realized, due to different types of light sources, various 
material textures, and the lack of a uniform material 
reflectivity of light. Tc greatly decrease display 
processing time, most of the image realism programs were 
implemented in hardware after procedure refinement. Thus, 
computer graphics realism has become a function of software 
time available or of the cost of complex hardware. 
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Figure 3 - WIBE-FBAME IMAGE 
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Figure u - UICDEN LIMES REMOVED 
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C. GRAPHICS SOFTWARE STRUCTURE 


In cider to provide a "state-of-the art", computer 
graphics software package, tie following capabilities were 
determined tc be the minimal requirements: 

1. Object to screen coordinate transformation; 

2. Clipping algorithm; 

3. Image scaling; 

U. Image translation; 

5. Imace rotation; 

6. Image shading; 

7. Hidden line removal; 

8. Hidden surface removal. 

These capabilities enabled the viewing of an image from any 
aspect and the generation of realistic displays in "real 
time". 


The data structure utuilized was selected for its ease 
of user implementation and its appicability to the input 
requirements of the eight procedures above. To explain the 
capabilities above (in Section III.), each pcint, or vertex, 
was defined as a one by four vector, such as: 

[x y z m], where m was a scale factor of the 3-D 
vector (normally m = 1) . 

This type of vector representation allowed all 
transformations to be defined as a four by four matrix. 
While the complete logic for matrix and vector utilization 
was provided in the Appendices of Ref. [1], a key advantage 
was that it facilitated the concatenation of 
transformations. Concatenation of matrices can be simply 
explained with the following example. 
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(here 1 =A B 


[x 1 y' z* 1] = [x y z 1] A 

[x» r 2- i] = [x* y z* i] b 

Can be equivalently stated as: 

[x» y» z» 1] = [ x y z 1 ] X, v 

In the next section, the algoritnms used to develops the 
graphics software package were briefly presented. since 
every display device has a smallest horizontal and vertical 
resolution size, the proper construction of an image 
reguired the computational utilization of these machine 
dependent features. An easy means to visualize this concept 
for any device was to construct a n by m dot matrix, where n 
was the number of horizontal lines and m was the number of 
dots per line, as shown in Figure 5. Thus, to display an 
image reguired that line segments be drawn tetween the dots. 
Similarly, it was necessary to determined the screens center 
in order to position the image in the middle of the display 
surface. The following terms were used consistently to 
define these guanties (see Figure 6) : 

V = m/2, was half of the horizontal resolution size; 
sx 

V = n/2, was half of the vertical resolution size; 
sy 

V = m/2, was the screen's horizontal center; 
cx 

V = n/2, was the screen's vertical center, 
cy 
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Figure 5 - N BY B DOT MATRIX SCREEN REPRESENTATICN 
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Eigure 6 - DISPLAY SURFACE DIHEHSIONS 
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III. GRAPHI CS SOFTWARE COMPONE NTS 


In this section the algorithm for each component was 
briefly presented along with any figures and general flow 
charts which aided with the explanation. For a detailed 
flow chart and Fortran program listing see Appendix A. 


A. OBJECT COORDINATE SYSTEM -THE DATA BASE 


The description of a 3-D object with this coordinate 
system allowed the user total flexibility in the selection 
of a convenient system of measurement and origin placement. 
Additionally, the user selected the viewpoint, which 
determined the viewing axis and thus, the initial viewing 
aspect. The viewing axis was the line defined by the 
viewpoint and the object system's origin. 

Image definition started with the specification of its 
vertices. As each vertex was input to the computer, it was 
assigned a consecutive index number. Similarly, as the 
other sets of image elements (edges, polygons, and polyhera) 
were input, they were consecutively indexed also. The x, y, 
and z object coordinate vales were stored in three real 
arrays, XE(i), YE (i), and ZE (i). Integer arrays were used 
to store indices decribing edges, polygons, and polyhedra. 
An edge was described by storing the index of one vextex in 
EDGE1 (i) and the second in EDGE2 (i) . A polygon was defined 
by storing the indices of the edges which composed its 
foundry in the array POLYGN (i, j) . The indices of these 
edges must be input so that: 
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{EDGE 1 (i) OB EDGE2 (i) } = {EDGE 1 (i +1) CB EDGE2 (i+ 1) } 

The polygons which described a polyhedron were input 
ccnsecutively to reduce storage requirements. Thus, a 
polyhedron was described by storing the index of the first 
and the last polygon in the array P0LYHE(i,2). This data 
structure allowed excellent imace flexibility, since any 
polyhedron could be rotated, scaled, or translated without 
altering the remaining display. 


E. IMAGE SCALING 


Ihe scaling transformation multiplied the object x, y, 

and z values by the scale factors S , S , and S 

x y z 

respectivley. If the scale factors used were not egual. 


image distortion resulted. 

[x« y 1 z* 1] = [x Y z 1] s. Where s 


3 0 0 

0 S 0 

y 

0 0 s 

z 

0 0 0 


C. TBANSLATION 


Ihe translation cf an image, in object coordinates, 
represents the physical movement of a 3-D object in the x, 

y, or z directions by the amounts T , T , or T , 

x y z 

respectively (or any combination of the three). The 
transformation T can be stated in matrix form as: 
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1=1 0 0 0 

C 1 0 0 

C 0 1 0 

-1 -T -T 1 

x y z 

D. IMAGE ROTATION 


The rotation of a 3-D object was broken into four 
distinct categories. Each involved the rotation cf the 
image through an angle, theta, about an axis in object 
coordinates. The following axes of rotation define the four 
categories: 

1■ Xr*B otatio n 


The image was rotated about the X-axis through the 
angle theta. The angle was measured in the clockwise 
direction abcut the origin, looking towards the origin from 
the positive X-infinity. The transformation, R , was 

defined as: 


0 cos 9 -sin 9 0 

0 sin 9 cos 9 0 

C 0 0 1 


2. IrAxis Rot atio n 


The image was rotated through the angle theta abcut 
the Y-axis, where the angle was measured as stated above in 
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transformation 


1. looking from positive Y-infinity. The 
matrix R , was defined as: 

y 


R 

y 


cos 6 0 sin 6 


0 1 0 
-sin 6 0 cos « 

0 0 0 


0 

0 

0 


3- irAxis Rot at ion 


The image was rotated about the Z-axis through an 
angle theta, which was measured as above from positive 

Z-infinity. The transformation matrix R was defined as: 

~z 

R = cos « -sin 00 0 

z 

sin 0 cos 00 0 

0 0 10 

0 0 0 1 

4. An Arbitrary Axis of Rot atio n 


The rotation of an image about an arbitrary axis 
required that the set cf transformations belcw be performed 
cn the image's object coordinates. Except for R^, these 

transformations were necessary to move the arbitrary 




into an axis for which a rotation algoritm (1., 2. and 3. 
above) already existed. 


The arbitrary axis was specified by any two distinct 
points cn it, (x,y,z) and (x',y',z'). The direction cosines 

for the axis, a , b , and c , were found by: 
x y z 
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a = (x'-x)/R, t = (y’-y)/R, c = <z'-z)/R, where: 
x y z 

R = V(x'-x)2 + <y'-y)2 + (z'-z)2 


Cne of the points, (x,y,z), was translated tc the 
origin using the transformation T, where: 

T = 1 0 0 0 

0 10 0 
0 0 10 
-x -y -z 1 


Now, the arbitrary axis was rotated into the Z-axis, 
by first rotating it about the X-axis through the angle 
alpha, as shewn in Figure 7. This transformation, R^, which 

placed the arbitrary axis in the X-Z plane, was defined as: 


0 cos Oi -sin <y 0 
0 sinCY cos O' 0 
0 0 0 1 


where: 


cosCY = c /v 
z 

sin(X - /v 

y 

v = Y't 2 +c2 


The arbitrary axis was then rotated about the S-axis 
through an angle beta, as shown in Figure 8, into the 
Z-axis. This transformation matrix, R , was described as : 
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-a 0 
C 0 


a 0" 

0 0 
v 0 


0 V- 


The rotation atout the arbitrary axis through the 
angle theta was now defined as a rotation about the Z-axis. 
Thus, the transformation, R^, was exactly as described in 3. 

above for E . The remaining task, was to return the 
z 

arbitrary axis, and the image, back to its original spacial 

position. This was accomplished by multiplying with the 

inverse of the transformations T, B , and E in the reverse 
~ “1 “2 

order. Using the principle of concatenation, the total 
transformation may be stated as : 


-1 -1 -1 
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Figure 7 - ROTATION INTO X-Z EIANE 
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Figure 8 - ROTATION INTO Z-AXIS 
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Figure S - EYE COORDINATE SYSTEM 
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Prior to this first coordinate transformation, the 
object coordinate values for the vertices were placed in the 
arrays XS(i), YS(i), and ZS (i). These arrays were used to 
generate the displayable information by the clipping, hidden 
line removal, or hidden surface removal algorithms. This 
prevented the original object's description from being 
destroyed or altered in these three procedures and allowed 
the graphics package to subsequently present different 
viewing aspects. Additionally, since clipping could remove 
an entire edge, the vertex indices were stored in EDGI (2,i) 
as: 

EDGE (1,i) = EDGE1 (i) and EDGE (2,i) = EDGE2 (i) 


The transformation matrices specifying V were formed 
as shown below when the viewpoint was located at (a,b,c) (in 
object coordinates) and the object was centered at the 
origin. Transformation T^ translated the viewpoint tc the 


origin by : 

1 = 10 0 0 
" 1 

0 10 0 
0 0 10 
-a -b -c 1 


A left-handed cartesian coordinate system was formed 

with T , where: 

“2 

I = fl 0 0 0 

"2 

0 0-10 
0 10 0 
C001 
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Figure 1C - VIEWING TRANSFORMATION 
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As shown in Figure 10, the system was rotated atout 

the y -axis through the angle theta which pointed the 
e 

Z -axis at the point (0,G,c). The transformation T was 
e "3 

specified as : 


T = ccs 0 0 

- 3 

0 1 

-sin 0 0 

0 0 

cos 6 = a/» , 
v = V a * + b 2 


sin 0 0 

0 C 

ccs 0 o 
0 1 


where: 


sin 0 = b/v 


Next, the coordinate system was rotated about the 

X -axis through the angle phi, as shown in Figure 11. This 
e 

pointed the Z -axis towards the object space origin, where: 


I 


u 


1 0 0 


0 , for: 


0 cos 0 -sin 0 0 

0 sin 0 cos 0 0 

_0 0 0 1 


cos 0 = v/(v z +c 2 ) , sin 0 = c/(v z +c z ) 


Ihis transformation arbitrarily selected the viewing 
axis as the line between the viewpoint and the object 

coordinate axis. It also placed the X -axis in the object 

e 

system's Z = c plane. Since the object coordinate system 

was user defined the logical initial viewing aspect was to 
look at the system's orgin cr center. The initial view due 
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to the position of the X -axis was acceptable and basically 
e 

irrelevant since the ixage could be rotated, translated, or 
scaled tc provide any desired viewing aspect. 
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Figure 11 - Z -AXIS HOTATED TOHABCS ORIGIN 
e 
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2. Eye to Screen C oordi nate Transformation 

This transformation completed the perspective 
projection cf the 3-E object onto the screen. As shewn in 
Figure 12, the display was generated by simply projecting an 
object's eye coordinates onto the plane cf the screeen.__ 

For a sguare display screen, one with equal horizontal and 

vertical resolution, the image was constructed without 

distortion, ty the following transformation: 

X =S (X/Z)V + V , where S = a/b 

s x e e sx cx x 

Y = S (Y /Z ) V + V , where S = a/b 

s y e e sy cy y 

Z = - 1/Z 
s e 

If the viewing screen was not square then S and S 

x y 

were modified as shown: 

If V > ¥ then S = (a/b) (V /V ) 

sx sy x sy sx 

If V > V then S = (a/b) (V /V ) 

sy sx y sx sy 

Ihis modification was necessary so that an un-distorted 
image cculd be displayed over an entire rectangular screen. 
Without it ere display dimension, typically the horizontal, 
would have teen elongated. 
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Figure 12 - PROJECTION OF DISPLAY POINTS ONTO SCREEN 
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DISPLAY CLIPPING 


The clipping procedure constructed a viewing pyramid 
which eliminated the undesirable effects of the perspective 
projection from object to screen coordinates, which were: 

1. points and thus objects may have been located behind 
the viewpoint; 

2. and objects may have exceeded the limits of the 
viewpoint (i.e. were located off the screen - were 
non-displayable). ^ 

The clipping of an image was performed on the image's data 
while it was expressed in eye coordinates to simplify the 
operation (as explained in Ref. [1]). 

As shown in Figure 13, the geometry cf the viewing 
pyramid dictated that for a point to be visible the 
following conditions must be satisfied: 

-z < s X S Z (1) 

e x e e 

-Z < S Y < Z (2) 

e ye e 

Thus, a transformation from eye coordinates to a 
"clipping" coordinate system was described as: 

[X Y 2 1 ] = [x Y Z 1] N , where: 

N = a/b 0 0 0 

0 a/b 0 0 

0 0 10 

0 0 0 1 
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Figure 14 - DISPLAX SCREEN DIVISION AND CODING 
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This coordinate system was established tc display only 
viewable points. Additionally, if an edge's endpcirt was 
located outside of the viewirg pyramid, this routine located 
a point on the edge which satisfied equations (1) and (2) 
above and became the new endpoint. By dividing the plane of 

the screen, the X -Y plane, into nine sectors (Figure 14) , 
s s 

the location of the twc verticies of an edge was determined 
using the inequalities (1) and (2). 

Since a vertex can be used as an endpcint of several 
edges, the clipping procedure placed the x, y, and z values 
of the two endpoints into the arrays X (2) , Y(2), and Z (2) , 
respectively, as each edge was examined. Thus, if new 
endpoints had to be computed for this edge, the values of 
the original verticies were not destroyed. Since Fortran IV 
did not support binary operations, an integer array, 
ICHK ( 2 , 4), was used to code the location of each endpoint as 
fellows: 

If X (i) < -Z(i) then ICHK(i,1) = 1 else ICHK(i,1) « 0 
If X (i) > Z(i) then ICHK(i,1) = 1 else ICHK(i,1) = 0 
If Y(i) < -Z (i) then ICHK(i,3) = 1 else ICHK(i,3) = 0 
If Y (i) > Z(i) then ICHK (i,4) = 1 else ICHK(i,4) = 0 

If both endpoints were displayable, no action was taken 
and the next edge was examined. If both vertices were to 
the right of the viewing pyramid (or both to the left, or 
botn above, cr both below), the entire edge was deleted from 
the display (i.e. if ICHK(1,j) = ICHK(2,j) = 1 then discard 
the edge). Until both endpoints were displayable or the 
edge could be rejected, new points had to be computed cn the 
edge. This computation has often been termed "pushing" the 
endpoint towards the display area. The pushing cf the 
endpoint was accomplished by utilizing the 3-D, parametric 
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representation of a line to select that point where: 


IX | 

c 


|Z | if inequality (1) was not satisfied; 
OH 



IZ | if eneguality (2) was not satisfied, 


If inequality (1) was violated, the following sets of 
equations are used to compute the new endpoint: 

If X (1) < —Z (1) : then: 

t= [ Z (1) +X (1) ]/ ([ X (1) -X (2) ]-[Z(2)-Z (1) ]} 

Z 11) = t * [ Z (2) -Z (1) ] + Z(1) 

X(1) = -Z(l) 

x (i) = t * [ x (2) -y (i) ] + y (1) 

If X(1) > Z(1) then: 

t = [Z(1)-X(1)]/{[X(2)-X(1) ]-[Z(2)-Z(1) ]} 

Z(1) = t * [Z(2)-Z(1) ] ♦ Z(1) 

X(1) = Z(1) 

Y (1) * t * [Y (2)-Y (1) ] + Y (1) 

When inequality (2) was violated, the equations used to 
compute the new point were those above with every X replaced 
with a Y and vice versus. 


G. HIDDEN 1INE REMOVAL 


The method utilized to remove hidden lines from 3-D 
objects was developed by John Warnock at the University of 
Utah. The program was interpreted from a SAIL program 
listed in Ref.£1]. This procedure required that the object 
coordinates be transformed to eye and then to screen 
coordinates without any intervening clipping of the image. 
The algoritm was broken into three main sections, the 
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looker, the Thinker, and the Controller. The storage cf the 
vertex indices in the array EDGE (2, j) Has re-structured so 
that the index of a polygon's first edge cculd be used to 
link to the index of its second, and its seccnd could link 
to the third edge, etc. Thus, an initializing subrcutine 
linked each polygon's edges in the array EDIINK(i). 

The concept of linked lists uses the index of edge i to 
produce the index of the next edge for the same polygon. 
Since an edge could be common to two polygon's, the first 
polygon to link the edge i found the storage location 
EDLINK(i) unused. This first polygon stored the index of 
its next edge, j, in ECLINK (i) . Because twc polygons can 
have at most one edge in common, the second polygon tc link 
edge i, could not use this same storage location, EDlIKK(i). 
Therefore, the number of edges, called edgem, was increased 
by one and the vertices for edge i were also stored as 
shown: 

EDGE<1,EEGEM) =ECGE(1,i) 

EDGE (2,EDGEM) = ECGE(2,i) 

Additionally, the storage location EC1INK(EDGEM) was 
used to stcre the index of the next edge for this second 
polygon. This ordering usually doubled the storage 
requirements for edge definition. The vertex indicies were 
further ordered in the array EDGE (i,j), so that: 

EDGE (2,j) = EDGE (1,j*1) 

The data was structured, using linked lists (integer 
arays) and pointers to the first element cf the list, as 
follows: 

1. POLLNK-a list of polygon indices ordered ty the 
polygon closest tc the viewer with pointer POLPTE; 

2. POLEDG- contained the indices of the first edge of 
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each polygon. The pointer used was the polygons index; 

3. EDIINK-contained the linked list of edges. The 
pointer to a polygons second edge was the index of its 
first (found in PCLEDG) ; 

4. POLLST-a list of polygons which were determined by 
the Locker to he either surrounders or intersectors 
(which are explained below). The pointers to the list 
were SUfiBND and INTER. 

An example of the usage of these linked lists was 
provided belcw; 

INTER was the last polygon (index) added tc the 
list of intersectors and thus, was at the head cf the 
list. 

POLLSI(INTER) = P, where P was the index cf the 
second polgon on the intersectors list. 

POIEDG(P) = E ^ , where E^ was the index cf the 

first edge for polygon P. 

EDIINK(E^) = , which was the index of the second 

edge fcr polygon P. 

EDGE (1 ,E ^) = and EDGE(2,E^) = V , where and 

were the indices of the two vertices discriting edge 


The concept of linked lists was utilized extensively for 
toth the hidden line and hidden surface removal routines. 
In this proceedure, a display window, which was initially 
the entire screen, was examined against each polygon. Each 
window could be classified as; 

1. nothing was contained in this window; 

2. the information contained in this window was simple 
and could be displayed; 
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3. or the information contained in the window was too 
complex. 

Situations 1. and 2. resulted in a successful processing of 
the window. The next window on the stack could then be 
examined. The last classification was a failure and caused 
the window to be divided into four windows of equal size. 
These new windows were then pushed onto the stack. 

The first important element needed to process a display 
window was the computation of the depth of the plane 
determined by a polygon at the four corners of the window. 
The planar eguation can be stated as: 

Ax + Ey + Cz ♦ D = 0 

The coefficients. A, B, C, and D, can be found from the x, 
y, and z values of any three pcints contained in the plane 
which are not colinear. Since the corners cf a window are 

specified as X and Y values in screen coordinates, the 
s s 

depth of the polygon was computed by simple substitution 
into the plane equation. 

The last important concept needed to process a window 
was the classification of each polygon as: 

1. an intersector of the window; 

2. a surrounder cf the window; 

3. or disjoint from the window. 

These concepts are clearly portrayed in Figure 15. The 
classification of all polygons was performed by the Looker. 

To determine whether a polygon was an intersectcr, it 
was sufficient to find any one of the polygon's edges which 
intersected the window. This determination was made by a 
clipping subroutine, which was very similar to that in F. 
If none of the edges intersected the window, the polygon was 
determined tc be a surrounder or disjoint from the window by 
computing the angle "about the window through which each 
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edge passed." ^ The sum of these angles, as all edges of a 
polygon were processed, vculd equal ±360 degrees, if the 

polygon was a surrounder of this window, as shown in figure 
16. The actual computation of each edge's angle was 
implemented ty dividing the X -Y plane intc nine sectors, 

as shown in Figure 17. The window was located in the center 

region, and the outer eight regions were numbered as 
shown. The endpoints were located, just as in the clipping 
routine in F. , and assigned the proper sector number. 

The edge's "angle" was the number of sectors which an 
edge entered, not counting the sector of the first endpoint. 
A polygon which surrounded the window had an "angle" of ± 8, 
and a disjoint polygon had a zero angle. Extremely complex 
polygons could have an angle equal tc ± 16, or higher 
multiples of eight, by surrounding the window two or more 
times. However, usage of such complex polygons was 
unnecessary to construct any image. Because the incremental 
angle ( ) was defined as the difference between the 

sector values of the two endpoints, one problem of computing 
an edge's angle occurred when the magnitude was greater than 
four. Since no linear edge could enter mere than four 
sectors (starting frem the sector of the first endpeint), 
the edge's angle was adjusted when the magnitude was greater 
than four by : 

If A 0( > 4 then A Of = AC* " 8 

if Ac* 1 < _i * then Aoc = A oi * 8 
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Figure 15 - POLYGONAL CLASSIFICATIONS 
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Figure 17 - DISPLAY AREA DIVISION AND CODING 
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The last angle computational problem cccured when the 
magnitude was egual to four, as shown in Figure 18. The 
correct sign in this situation was not rendered by taking 
the difference of the two sector values. Counter-clockwise 
rotation or movement about a polygon's edges should have 
yeilded a positive result. The problem was resolved by 
selecting ary point between the two endpoints which was not 
in either of the vertices' sectors, as shown in Figure 19. 
Ey dividing the edge at this point, the correct angle could 
be computed ty summing the angles of these "two edges". 
This angular computation for an edge was determined in the 
clipping subroutine used by this proceedure. 

The linked lists of the classified polygons were then 
passed to the Thinker. The surrounder list was processed 
first to determine which pclygcn was closest to the viewer 
by computing the depth of each of these polygons at the four 
window corners, as shown in Figure 20. Provided the closest 
polygon, called the hider, was not penetrated by ar.cther 
polygon, these four depths were used to determine if an 
intersector polygon was located completely ir front of the 
hider within the confines of the display window. If an 
intersector was completely hidden from the viewer by the 
bider, it was removed from the list. If the final 
intersector list contained only one polygon, then that part 
of the polygon's edges which were inside the window were 
displayed. If the list contained more than cne intersector, 
or if any intersectcr polygon penetrated the plane cf the 
hider, the Thinker announced failure for that window. if 
the hider was penetrated by another surrounder, the Thinker 
announced failure before examining the intersector list, 
lhe penetration cf cne polygon by another was shewn in 
Figure 21. Whether polygon B would be classified as an 
intersectcr or a surrounder penetrating polygon A would 
depend on the placement of the window. When the complexity 
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of the display could not fce resolved and the size cf the 
window had been reduced to the display device's smallest 
resolution, a dot was displayed at the window's lower, left 
corner. In this manner, the penetration of cne polygcn by 
another, which described a line, was displayed as an implied 
edge. 


If failure was announced, and the size cf the window was 
larger than the smallest resolution, the display window was 
divided into four equal windows and pushed onto the tcp of 
the stack. The Controller then selected the next window on 
the stack and passed it to the Lcoker. This entire process 
was summarized by the flow chart in Figure 22. 
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Figure 18 - ANC-LE MAGNITUDE Cl FOUR 
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Figure 19 - EDGE DIVISION FCR PROPER ANGIE CONFUTATION 
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Figure 20 - 



DEPTH COHPUTATION CF POLYGON AT FOUR CORNERS OF 
WINDOW 
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Figure 21 - POLYGONAL PENETRATION 
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Figure 22 - HIDDEK LINE REMOVAL FLCH CHART 
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B. HIDDEN SURFACE REMOVAL 


The algorithm utilized for this procedure was developed 
by G. S. Batkins at the University of Utah. Ihe program was 
interpreted from a SAIL program presented in Ref. £1], 
Nhile the hidden lice removal algorithm concentrated on 
linked lists of polygons, the hidden surface algorithm 
processed the display with lists of edges. Additionally, 
the vertices had to be expressed in clipped, screen 
coordinates. Although the user still generated the data for 
an image as stated in A., this procedure displayed the 
polygonal surfaces as a solid plane using shading or colors. 
Ihus, a 3-D object should become much more realistic when 
displayed with shaded surfaces vice wire frames. While the 
algorithms presented previously in section III., can be 
utilized on any display device, this concept was developed 
specifically for raster scan CRT's. 

A raster scan is a special type of CRT, which is very 
similar to the television in most homes. The vectored CRT's 
and the direct view storage tunes generate a display by 
pointing an electron beam to a desired location on the 
display screen and then moving it to any other screen 
location. This process illuminates the phosphorous screen 
to produce a single line segment. The typical television 
receives an analogue broadcast signal which generates a 
single ^horizontal line of the screen's image at a time. At 
the end of each line a horizontal sync pulse is received to 
move the electron beam down one line, and to the left-hand 
edge. When the last horizontal line has been displayed a 
vertical sync pulse moves the beam to the top, left-hand 
corner of the screen. Since the phosphorous screen remains 
illuminated a very short time, the image must be constantly 
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refreshed, typically at a rate of thirty times per second. 

A raster scan display device receives its image (and 
refresh) information from random access memory (RAM) refresh 
planes where the image is stored as a sequence of individual 
tits. Each bit of a memory plane determines the 
illumination of a single element on one horizontal dislay 
line (also called a scan line or a raster). A picture 
element, called a pixel, is the smallest screen resolution 
size. The standard sixteen level grey shading requires four 
bits, one bit on four planes, to represent the shading of 
one pixel. Similar memory requirements are needed to 
display an image with sixteen colors. While the vectored 
CRT's have acheived resolutions on a display screen of 4096 
lines with 4096 elements per line, the finest resolution 
available with raster scan devices is 1024 by 1024. Thus, a 
sixteen color, raster scan display with high resolution 
required fcur million bits of HAM. Because of this 
extensive memory requirement, the developement of this type 
of display device followed that of the small, lower cost 
electronic memory. 

As shewn in Figure 23, the intersection cf the plane of 
a scan line with a polygon was a line segment. (Scan line k 
corresponds to the Y = k plane.) This line segment's 

endpoints were defined by its Xleft, Zleft, Xright, and 

Zright values, which were the X and Z coordinates of the 
s s 

intesecticn of the scan line with two cf the polygon's 

edges. The two, 2-D equations below were used to find the 
intersection of an edge with each scan line. 

X = a Y + b (1) 

s s 

Z = c + d (2) 
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Figure 23 - SCAN FINE INTERSECTION OF A POLYGON 
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Figure 24 - SEGMENT CLASSIFICATION 
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The coefficients, a, b, c, and d, were guickly obtained 
using the eguations for a 2-D line with the coordinates of 
the two vertices describing each edge. 

With these eguations the X and Z values of the 
intersection of an edge with scan line k+1 was the values at 
scan line k plus their respective slopes, a and c. Since 
each edge, and thus each segment, sustained only an 
incremental change between scans lines, the display was also 
assumed to remain constant between scan lines. This scan 
line coherence of the display was used to decrease the time 
reguired to process an image. 

To process each scan line, it was necessary to divide 
the line intc spans which could be more easily resclved. 
The content of a span cculd be categorized as shewn in 
Eigure 24 and as described below: 

1. "The span contained only one segment. 

2. One segment was closer to the viewer than all 
others and it was a spanner. A spanner was a segment 
where Xleft < Span left and Xright > Span right, as 
shown in Figure 24. 

3. There was a simple intersection of the only two 
segments in the span, and both were spanners. This 
span was divided at the intersection intc two spans and 
processed as in 1. 

4. The display was too complicated in this span so it 
was divided at the left-nost segment endpoint, or at 
the spars mid-point if there was no endpoint. The new 
spans were then processed." ^ 

Since the lower, left-hand corner of the raster scan 
machine at the Naval Postgraduate School, the RAMTEK, was 
indexed as (0,0) each edge was ordered so that the index of 
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the vertex with the largest Y value was stored in 

EDGE (1,i) . The integer value formed by truncating this 

vertice's Y value determined the first scan line that an 
£ 

edge would enter the display. A linked list cf the incicies 

of the edges which ertered on each scan line was stored in 
the array ENlLST(i). The index of the first edge to enter 
cn scan line k was stored in YENTER(k). 

As each edge entered the display, the X and Z values and 
the scan line coherence factors, the slopes for equations 
(1) and (2), were computed. Since the object of this 
algorithm was to display polygonal surfaces, the indices of 
the current segments of a polygon were linked in POISEG (i). 
The segments were ordered by increasing X values of their 
left endpoint and the first segment's index was placed in 
SEGLST(p) for polygon p. A segment's index pointed to a 
block of storage which defined the endpoints X and Z values 
in the arrays XL EFT (i) , ZLEFT(i), XRIGHT(i), and ZRIGHT (i) 
and their respective slopes DXLEFT(i), DZLEF1 (i), DXRGHT (i), 
and EZRGHT(i). Integer arrays IYLEFT (i) and IYRGHT (i) were 

used to indicate when an edge, the source cf one endpoint of 

a segment, was exiting the display. 

To properly insert an entering edge into a polygon's 
segment lists required the comparison of the edge's X value 
for this scan line to the Xleft and Xright values of all cf 
the active segments. If two edges entered cn the same scan 

line at the same X coordinate, as shown in Figure 25, the 
s 

edge with the largest slcpe was entered first. Thus, to 

enter the twc edges on scan line k+1, between the two 

existing edges, a new block of storage was added for each 

entering edge. As each edge exited the scene, it was 
removed from its half of the storage block, as shown on scan 
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line k+2 in figure 25. After all additions and deletions 
had been performed, the list of segments was sorted to 
conslidate storage. This entire process has been portrayed 
in Figure 26. As stated before, an edge usually separated 
two polygons. To eliminate redundant operations, the array 
E(2,i) was used to store the indices of the polygons common 
to edge i. fchen an edge entered the display and separated 
two polygons, its values of intersection were added to the 
blocks of storage for both segment lists. 

This algorithm was divided into the same three parts, 
the Thinker, the Looker, and the Controller, as the hidden 
line routine. The Looker compared all segments which 
intersected a span and developed sufficient information for 
the Thinker to process it. Provided the contents of a span 
satisfied the categories 1 through 3 above, the Thinker was 
able to generate the data required for displaying this span. 
If the information contained in the span was too 
complicated, it was divided by the Controller. The 

sucessful scan line division points normally occurred at the 
left-most endpoint of a segment. since this division 

point's location on the next scan line can be predicted, it 
was stored and used to decrease the time required to process 
the entire image. The computation of the X and Z values of 
a segments endpoints and an updated, sorted list of segments 
was also performed by the controller. This active list of 
segments was sorted in the Xsort lists, IXSLFT(i) and 
ZXSBGT(i), by increasing X values. These twc lists provided 
the index of the segment to the left and to the right of 
segment i by: 

k = IXSLFT(i) : was the index of the segment to tae 

left. 

j = IXSBGT(i) : was the index of the segment tc the 

right. 
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Figure 25 - ENTERING AND DEPARTING SEGMENTS 
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Figure 26 - UPDATE OF SEGMENT BLOCK STORAGE 
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The last bookkeeping task performed by the Controller 
was to divide the list of segments into the following 
categories: 

1. "SEGOUT- the right edge of the segment was 
contained in this spap. (This segment did not appear 
in any span to the right of this one.) 

2. SEGACT-the right edge of the segment extended 
beyond the right limit of this span." ^ 

If this span was displayed, then the SEGCDT list did not 
need to be considered in subseguent spans and could be 
discarded (until the next scan line). The SEGACT list was 
automatically added to the next spans active segment lists. 
If the span failed to be processed, the two lists were 
combined and compared to the new span. 

When the segment lists were passed to the Looker, the X 
and Z coordinates of the left-most and the right-most parts 
of a segment in this span were computed (see Figure 27 where 
the following terms are portrayed: sxleft, szleft, sxright, 
and szright). A box was constructed about the first segment 

examined in the X -Z plane which entirely surrounded that 
s s 

part of the segment which intersected the span. As each new 

segment was compared to the box, the box was enlarged to 
include it. If a segment completely hid the box from the 
viewpoint or if it was a spanner, then the count of the 
segments in the box was reset to one and the tox was made to 
enclose only this segment. This was shown in Figure 28 
along with the definitions of the box X and Z limits, 
txleft, bzleft, bxright, and bzright. 
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Figure 27 - SPAN DESCRIPTIVE TERMS FCR A SEGMENT 
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Figure 28 - SEGMENT BOX DEFINITION 























The information passed to the Thinker was box count and 
box type. If the box count was zero, the Thinker did 
nothing and the Controller began processing the next span. 
A box count cf one indicated a single segment existed in the 
span. The segment*s sxleft, sxright and index were stored 
with the results of previous spans so that the entire scan 
line could be displayed at one time. Hhen the box count was 
greater than one and the box type was equal to one, the span 
contained a simple intersection of two spanners. Both 
segments' sxleft and sxright values and their indices were 
stored with those of previous spans for this scan line. If 
the box ccunt exceeded one and the box type equaled zero, 
the Thinker announced failure. The Controller then divided 
the span and began to process the left half cf the old span. 

An important element of all hidden surface or line 
elimination algorithms has been the proper display cf an 
implied edge. Hhen the intersection of two polygons caused 
an implied edge, this proceedure stored the value DIV, the 
division point caused by two intersecting spanners, cn its 
first occurence (see Figure 29) . The occurrence of this 
division on a second scan line provided sufficient 
information to compute its X and Z scan line coherence 
factors. The implied edge was then added to a dummy segment 
block and sorted with the other segments in the Xsort lists. 
However, this dummy segment block was net passed to the 
looker, but was used to divide the scan line into the proper 
spans. When the division point due to an implied edge 
divided the same segment, the implied edge was discarded. A 
generalized flow chart was provided for this algorithm in 
Figure 30 to summarize this entire section. 
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Figure 29 - IMPLIED EDGE GENEB&1I0N 
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figure 30 - HIDDEN SURFACE REMOVAL fLOW CHART 
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I. IMAGE SHAEING 


An extremely important aspect of computer image realism 
was the generation of an appropriate shading algorithm. The 
realistic algorithms required complex software and thus, 
uore computation time. Another alternative used has been to 
iaplemenx the algorithm with sophisticated and expensive 
hardware. Although this was one key element to a 3-D 
graphics language, a software implementation would have 
required half again as much reasearch and time. Therefore, 
this aspect was left for future developement. The usage of 
colors or shading for object definition was input via the 
subroutine 1NITAL along with the ether image data (a 
complete description was included in Appendix A). 
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IV. HABDWAEE AND SO FTW ARE CONSIDERAT IONS 


The realistic implementation of a standard graphics 
language throughout a large crganizaticn would require that 
the majority of the software could be input to all host 
computers with minimal alterations required for specific 
display devices. Therefore, the primary intent was to 
reduce display device dependence to the fewest number of 
subroutines. The processing of a 3-D display in real time 
was another important consideration. In the following 
sections the algorithms in III. were divided intc four 
groups and their implementit ion on the display devices was 
presented. 


A. IMAGE DISPLAYED - ALL LINES SHOWN 


The device dependent software which had to be utilized 
to display an image and the different screen cr display area 
of each machine were the two non-portatle aspects of this 
entire graphics package. Since the intent of this effort 
was to produce a language which appeared to be device 
independent to the user, the problem of various display 
areas was resolved in the image data input subroutine, 
INITAL. When more than one display device was supported by 
a host computer, the user had to select the appropriate 
display device number. INITAL then chose the correct line 
and element resolution and the location of the display 
area's center. The four devices used in this project were: 

1. 1EK1R0NIX 4012- a direct view storage tube; 
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had 


2. ADAGE AGT-10- a vectored CRT; 

3. VEHSATEC- a hard copy device which also 
electrostatic shading capability; 

4. BAH1EK- a raster scan CRT. 

The format of the data input by the user was stated in 
the graphics language description preceeding INITAL. The 
required data included the image description as stated in 
III. A. and the following information: 

1. a - the distance from the screen tc the viewer; 

2. t - the vertical dimension of the screen; 

3. the x, y, and z coordinates of the viewpoint; 

4. the select display device number; 

5. the index of the color table to be utilized (when 
using the RANTER); 

6. and the index of the color or shading for each 
polygon. 

A complete explanation of the input data was provided in 
Appendix A with the program listing of INITAI. 

To display the image as input without removing hidden 
lines, the user then had to call the subroutine DISFLY. 
EISPLY called the subroutines listed below. 

1. RDYC1P- performed the object to clipping coordinate 
transformation on the vertices and stored the results 
in the arrays XS (i), YS (i) , and ZS(i); 

2. CUE- clipped the image against the viewable 
display area; 

3. SCBN- converted the clipping coordinates to screen 
coordinates. 

The image could now be displayed by the selected device. 
First, each machine had to be initialized by a single device 
dependent subroutine call. Next, the line segments of each 

86 



edge, also called vectors, were drawn by another device 
dependent subroutine. Finally, all devices, except the 
BAMTEK, reguired a subroutine call which terminated the 
image. The ADAGE also reguired an image subroutine which 
developed a display list of vectors which was used to 
refresh the display screen. Thus, the device dependent 
portions for direct image viewing could typically be 
contained ir three subroutines and at most four when a 
display list was reguired. The display device was 
initialled by the subroutine INITAL and vectors, or line 
segments, were drawn by DISPLY. The graphics software 
package was always terminated by a call tc FINISH, which 
when required called the device dependent subroutine to 
terminate the display. 

The only calling parameter which was reguired for EISPLY 
was a two element integer array, IE (2), which was used tc: 

1. display a single polyhedron by setting 
IR (1) = IB(2) = the index of the polyhedron; 

2. display a group of consecutively input polyhera, 
where: 

IB (1) = the index of the first; 

and IR(2)=the index of the last pclyheron; 

3. or display the entire image as input in INITAL by 
setting IR (1) = 31. 

Except fcr the TEKTRONIX, the display cf a wire frame 
image by video graphic machines was performed in "real 
time". Beal time was Defined as a period cf time so short 
that the user could not detect a time lag between the 

programs execution and the complete drawing of the image, 

iihile the computation of the image require an insignificant 
amount of time for the host computer of the TEKTRONIX (an 

IBM-360), the vector display speed of this device was very 

slow. The V EBSATEC was also a relatively slew device, but 
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hard copy machines have not been expected tc produce "real 
time" displays. 

The construction of this general graphics language did 
emit some specific device versatilities. The VERSATEC had 
the capabilities to re-define its display area, resolution 
size, and line thickness, and could produce shaded images. 
These capabilities could have been added, but were left for 
subsequent efforts due to time constraints. 


E. INTERACTIVE SOETHARE AND HARDWARE 


The transformations for rotation, scaling, and 
translation provided the means fer complete image movement 
in three dimensions. These procedures, which have been 
implemented with hardware at some installations, coupled 
with interrupt devices, such as alphanumeric keyboards, 
function switches, graphic tablets, joysticks, light pens, 
and track balls, provided a user with a complete interactive 
viewing capability. The subroutines written to rotate, 
scale, and translate a single pelyhedren or the entire image 
were completely device independent. 

The calling parameter IR(2), which was used by each of 
the three subroutines, was defined and utiilized as stated 
in A. above. Thus, new viewing aspects could be generated 
for a single polyhedron, a group of objects, or the entire 
displayable image. The usage of the scaling subroutine also 
required that the user provide the scale factors for the x, 
y, and z coordinates. The additional calling parameters 
required for image translation were the distances in the x, 
y, and z directions which the 3-D object was to be moved. 
The sign required for these distances was opposite to that 
of the standard velocity vector describing the objects 



motion in this direction. Rotation calling parameters 
included the number specifying the axis and the angle of 
rotation. If an arbitrary axis was selected, the 
coordinates cf two distinct points had to be passed into the 
subroutine also. The usage of these routines was defined in 
the comment section preceeding INITAL. The three software 
transformations, rotation, translation, and scaling, 
reguired a minute amount of computation time. Any 
combination of a few of these three called between INITAL 
and DISP1I cculd be performed in "real time". 

Direct view storage tube display devices have an 
extremely limited interactive capability due to the method 
used to clear the screen. This device had a writing cathode 
which traced the image on a fine wire mesh which was lccated 
just behind the phosphorous screen. Initially, the entire 
wire mesh was negatively charged. Vectors drawn on the mesh 
by the writing cathode caused these line segments to become 
positively charged. These positively charged areas 
accelerated and passed the electrons emitted from a second 
cathode, which was issuing a "flood" of electrons to refresh 
the image or. the phosphorous screen. To clear an image, a 
large positive pulse was applied to the wire mesh. This 
caused a large flash to spread across the screen. Since the 
flash disrupted any possible display for several seconds, 
the usage of this type of device for rapidly changing, 
interactive images was highly unrealistic. 

Vectored CRT's have acheived an extremely high degree of 
resolution and support most interactive interrupt devices. 
With this type of display machine, the image on the 
phosphorous screen was refreshed by storing the entire image 
in a vector list. A highly complex, static image, one 
containing several thousand vectors, could cause the display 
to begin flickering. This type of CRT can produce only a 
limited number of vectors before the phosphorous illuminated 
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for the first vector begins to dim. Thus, shaded images, 
which would require many vectors for even a simple 3-D 
object, can not feasibly be produced using a vectored CRT. 
However, excellant 3-D graphs and extremely complex 
wire-framed objects have been visulized, using multiple 
colors, with these devices. 

The ADAGE AGT-10 was an extremely versatile device and 
had an alphanumeric keyboard, function switches, function 
knobs, a joystick, and a light pen as interactive 
capabilities. Since this machine was operated in a 
stand-alone mode, these interactive interrupt devices were 
easily utilized through the users application program, 
except the light pen. The light pen could only be accessed 
and utilized in the image subroutine. This extremely, 
device dependent capability was not included in the graphics 
software, since its usage would have required a through 
knowledge of this machines software for even a simple 
application. 

The BAHTEK's interactive devices consisted of an 
alphanumeric keyboard and a set of function switches which 
could be used to position a cursor. The cursor's screen 
coordinates were obtainable through device related software. 
The EAMTEK's host computer was a PDP-11, which primarily 
supported the software "C". Since the graphics software for 
the RAHTEK was written in C and there was no software 
interface written for Fortran IV, only those subroutines 
required to perform hidden surface elimination were 
translated to Fortran. Therefore, the interactive 
capability of this device was limited to input via a host 
computer's alphanumeric terminal. The resolution of this 
video graphics device was 240 lines with 640 elements per 
line. Its lack of vertical discrimination provided poor 
image continuity in this dimension. Since this was the only 
raster scan device available at this school, its utilization 
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was necessary to implement the hidden surface algorithm, 


C. HIDDEN LINE REMOVAL 


The hidden line removal algorithm described in Section 
III. was invoked by calling the subroutine REMOVE after 
INIIAL and any desired interactive subroutines. REMOVE 
called the subroutines RDYCLF, SCRN, WARNCK, and DISEI2. 
Ihe calling parameter required by REMOVE was the integer 
array IR(2) r which was described in A. Since this algorithm 
clipped each edge against display windows to process an 
image, the vertices of the pclyhderon to be displayed were 
passed tc WARNCK expressed in un-clipped screen coordinates. 
WARNCK contained the Looker, the Thinker, and the Controller 

described in III. As the display was processed, the X and 

s 

Y coordinates of the two endpoints for each vector were 
s 

stored in two arrays. This storage reduced the number of 

device dependent subroutines added to the graphics package 
by this algorithm to cne, DISPL2. DISPL2 generated vectors 
for the display exactly as performed by DISPLY. The image 
subroutine used by DISPLY to created the vector list fcr the 
ADAGE was also used by this subroutine. 

The algorithm as presented in Section III. displayed 
each edge of a simple polyhedron, like a cube, by failure. 
Display by failure means that the Thinker was unable to 
resolve any display window and a dot was displayed when the 
window's size was reduced to the smallest screen resolution. 
An edge was displayed as a line cf dots and was extremely 
grainy. This consistent failure had occurred because each 
edge was common to two polygons. When the display window 
was reduced so that it contained only cne edge, the 


91 



intersectcr list still contained the indices of two 
polygons. Ihe Thinker announced failure and the Controller 
divided the window. The computation time required to 
process and display a cute (which has a maximum of three 
viewable surfaces and nine viewable edges) exceeded twelve 
minutes of CPU time on the IBH-360. Additionally, almost 
9,000 storage locations were required to store the endpoints 
cf these single dot vectors. 

To reduce the occurrence of display by failure, the 
intersector list was not rejected if ctly two polygons 
remained after comparison with the hider. When the list 
contained two indices, the number of edges which intersected 
the current display window was determined for both 
polygon's. Provided there was only one edge for both, the 
vertex indices of one polygon's edge were compared to those 
of the seconds to ensure it was the same edge. When this 
procedure found a common edge, the endpoints cf the vector 
intersecting the window were stored. This addition to the 
Thinker reduced the CPU time to approximately 20 seconds and 
the storage requirements to about 300 locations. 

Even though the display time had been reduced by a 
factor of thirty, this hidden line removal procedure was not 
even remotely acceptable for a real time display. While a 
more complex Thinker could possibly reduce the computation 
time, hardware implementation remains the only feasible 
method for a real time display using this algcritm. 
Additionally, the realism of wire-frame images, even with 
hidden lines removed, was marginal at best. As stated 
previously, shaded surfaces, or solid images, can net be 
displayed with this type of CE 3. Therefore, realitic 3-D, 
"real time" computer graphics must be performed using a 
raster scan CBT with a hidden surface removal algorithm. 
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HIDDEN SURFACE REMOVAL 


The hidden surface algorithm was implemented by calling 
the subroutine SURFAC after INITAL and any of the desired 
interactive procedures. As with all ether procedures in 
this package, the hidden surfaces may be removed and the 
display created for one or a set of polyhedra or the entire 
image input to INITAL as determined by the two element 
integer array IR (2). The object coordinate data was 
transformed to clipped, screen coordinates by calling the 
subroutines RDYCLP, CLIP, SCRN, and SHOWIN. The remaining 
portion of this subroutine, SURFAC, was tbe Controller as 
described in Section III. The Looker and the Thinker were 
contained in the obviously named subroutines LOOKER and 
THINKER. 

The hidden surface removal algorithm added two device 
dependant subroutines to the graphics software package, 
which were: 

1. SHOSIN - constructed the desired color table and 
corrected the vertical coordinates; 

2. SHO k - displayed each scan as its image was 
resolved. 

The hidden line subroutines stored all of the image vectors 
until the entire display was resolved. However, even a 
single polyhedron would generate such a large list of 
vectors (or segments) to display solid surfaces that the 
storage requirements of a moderately complex scene wculd 
exceed realistic limits. 

The RAM1EK, and most raster scan CRT's, could produce 
sixteen intensity levels fer each of tfce three primary 
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colors. Thus, with all possible combinations of the shads 
of red, green, and blue, this device was able to display any 
16 of the 4096 colors at one time. Since a shading 
algorithm was not implemnted, SHOWIN was used to construct a 
color tatle, containing sixteen colors, which was needed for 
the applications program presented in Section V. 
Additionally, the RAMTEK's vertical scan lines (240 lines 
total) were twice as wide as each horizontal element of 
resolution. To prevent this rectangular picture element 
from causing image elongation, the vertical dimensions, V 

sy 

and V , were doubled. This procedure caused the 
cy 

horizontal, X , coordinates to be clipped with Y values 
s s 

which were twice their actual size. Thus, the large 

disparity between vertical and horizcntal resolution 
dimensions was not allowed to cause image elongation. 

SHOWIN was also used to divide all X coordinates ty two 

s 

before the controller began processing scan lines. Had 

these coordinate values not been divided, each scar line 
displayed would have been resolved and written intc BAM 
twice. 

The hardware used to input the image data onto the 
memory planes by the RAMTEK was called a "vector generator". 
Although the "vector generators" used by current raster scan 
CRT's disply vectors at speeds only limited by the memory 
plane write times, the RAMTEK's generation cf line segments 
was noticeably slow. Except for this slew vector display 
capability, this software algorithm developed shaded 
surfaces in "real time". Provided the user defined the 
shading or colors of each polygon, this type of display 
algorithm was shown to be a valuable tcol for rapidly 
changing, realistic presentations. 
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V. A THREE-L IHE SN ION AL GRAPHICS APP LICATION 


In order to demonstrate the capabilities of the 3-D 
Graphics Package an application program was written, which 
was motivated by the display of the torpedo test area at 
Keyport. Each test area, which was irregularly shaped, 
could be described by one or more convex polyhedra. A 
polyhedron's upper and lower surfaces represented the Puget 
Sound's air/water bcundry and its mud bottcx, respectively. 
Rotation, scaling, and translation of these polyhedra 
provided any desired viewing aspect. Since the display 
device at Keyport, the GENISCO GCT-3000, was a high 
resolution (1024 by 1024) display device, the hidden surface 
algorith was selected to present a realistic 3-D display. 
However, the Range Safety Officer needed tc see insice the 
polyhedral approximation of the torpedo test area not the 
closest polygonal surfaces. In view of this common type of 
display requirement, the algorithm was altered so that the 
hidden surfaces were displayed and the closest surfaces were 
deleted. 

A real time display of torpedoes and other test vehicles 
in the range was required to prevent them from running 
aground or leaving the test area. Hence, an algorithm was 
developed which detected the polygonal surface penetration 
by a vehicle. Obviously, detection of an accident after 
cccurrence wculd be absurd. Thus, an actual implmentation 
of this type of display would require the construction of a 
safety factor by displaying surfaces which were withir the 
actual test area boundries. This construction would be 
performed by the definition of the topedo test area when 
input to INITAL and actually did not alter the display 
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requirements. 


The location of a vehicle in the test area was provided 
by acoustic line of bearings which were used to determine 
its position. For this application, the fix provided by the 
bearings was assumed to generate the x, y, and z object 
coordinates required for this graphics package, since the 
acoustic information provided about the vehicle was not 3-D, 
the torpedo’s location was defined as a single point. 
However, the actual display of the vehicle was 3-D which 
provided the correct relative direction of motion 
information (called target angle) to the observer. Although 
all surfaces of the tcrpedo were the same color, its hidden 
surfaces were removed for display since this algorithm 
required such minute computational processing time. 
Einally, tne torpedo’s track, its last five positions, was 
displayed as a line in the vehicles color. 

A. FLAN AG SURFACE PEN ETHATION 

The theory used to determine if a point had penetrated 
one of the test area's boundries was originally conceived by 
I. G. Roberts to remove hidden lines from 3-E figures. it 
used the coefficients of the plane equation presented in 
Section III. G. "These coefficients, in the vector form [a 
b c d], were also the expression for a hciogeneous vector 
normal to the plane (homogeneous refered to the 
representation of a 3-D pcint as a 1 by 4 vector, where d 
was an arbitrary scale factor). If the dot product of this 
normal and a vector in the viewing direction was positive, 
then the polygon was a back face of the polyhedron and thus, 
non-viewable. Face normals were also used to compute 
shading paraneters. " 
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Additionally, if two points were on the same side cf the 
plane of a polygon, the dot product of either point with the 
normal vector would have the same sign. Since concave 
polyhedra would allow two interior points to be on different 
sides of a plane determined by a polygon, only convex 
polyhedra were utilized. Two adjacent ccnvex polyhedra 
could have at most one common polygonal face. Since 
penetration cf this common surface by a vehicle would have 
falsely indicated danger, these faces were eliminated. If a 
vehicle penetrated a test area boundry, its color was 
changed to red. 

E. HIDDEN SUBFACES DISPLArED 


To present the interior of a 3-D object, it was 
necessary tc eliminate the normally viewable surfaces. The 
hidden surface algorithm, which originally eliminated back 
planes, was easily modified to display the hidden surfaces. 
Only the Looker subroutine had to be altered. The Looker in 
Section III. compared all segments which intersected a span 
in order to locate one segment which was closest tc the 
viewer and hid all other segments in the span. This 
subroutine was altered to search for the farthest segment 
from the viewer which hid all other segments from a 
viewpoint located on the other side of the crigin (on the 
same viewing axis), as shown in Figure 31. 

Since the majority of this image was static, the tcrpedo 
test area cnly reguired processing to display the hidden 
surfaces initially and when the viewing aspect was changed. 
However, any moving image presented a special problem on 
this type of display device. To project the concept of 
motion on a raster scan CRT, it was not sufficient to simply 
translate or rotate and then display the new image, because 


97 



the original figure was written in the memory planes and 
would be displayed until replaced. With the GENISCO, it is 
possible to sample the memory planes defining the color of 
any one pixel, but the P.AMTEK, which was designed and built 
ten years before the GCT-3000, did not have that capability. 
1c delete an image of a torpedo with the GENISCO, one method 
would be to sample the two colors displayed cn the screen at 
the two endpoints of each line segment which defined this 
figure. Provided colors at these two points were the same 
and the vehicles image was small, this colcr could be stored 
for this segment. It could then be used to restore the 
original background when the torpedo's position changed. 
However, the time required for these background color 
computations could be more than that to re-process and 
display the entire image again. without this option, using 
the RAMTEK, the static image was input to the hidden surface 
display algorithm and re-drawn each time the torpedo's 
location changed. The torpedo's 3-D form and its track 
where then displayed at their new screen locations. 

To provide the maximum display flexibility, it was 
necessary tc provide a means to both relieve and display 
hidden surfaces. Thus, an additional calling paraireter, 
ILOOK, was added to the subroutine SURFAC (the Controller). 
Ihis parameter was used to enable the Controller to select 
the correct looker subroutine, where: 

1. ILOCK = 1, called LOOKER which removed hidden 
surfaces; 

2. and ILOOK =2, called LOOKR 1 which displaced the 
back pclygonal planes. 
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Figure 31 - VIEWPOINT REQUIRED TO DISPLAY THE BACK SURFACES 
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TOBPEDC JEST ABEA SIHULATION 


The torpedo test area utilized for this application 
program was constructed with three, box type polyhedra. The 
two end boxes had five polygonal surfaces and the middle had 
four. The test area, which was shown in Eigure 32 as a 
wire-frame image, had all surfaces displayed as light blue, 
except the bottom which was light brown. The torpedoes were 
colored black. 

Interactive input to this program provided new torpedo 
positions in x, y, and z object coordinates, which simulated 
periodic acoustic fixes and subsequent display update. Any 
two, distinct vehicle locations described its direction 
vector. Kher a new position was input, a new direction 
vector was formed between this and the last torpedo 
location. Tc compute the correct target aspect, these two 
direction vectors were used to find the angle through which 
the vehicle tad rotated, as shown below: 

cos « = (J b) / (||a|| ||b||) , 

where a and b were the two direction vectors and the 
||a|| operation represented the magnitude of vector a. 

The vehicle's image was then rotated about the axis 
which was rormal to both vectors. This normal vector was 
found by forming the cross product of the two direction 
vectors. The two points used to specify this arbitrary axis 
required by the subroutine ROTATE were the vehicles present 
position, (x y z 1), and a point found using the normal 
vector. These two operations were summarized algebraically 
as : 
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N = a a b , where g represents the operation of 
vector cross product. 


and 

<x* y* z* 1) = (x y z 1) + N 

Thus, the correct vehicle aspect was obtained by 
translation cf the image to (x y z 1) and rctation about the 
arbitrary axis specified above through the angle theta. 

lo realistically describe a 3-D object by planar 
polygonal mosaics was difficult if the object was composed 
of curved surfaces like a torpedo. Its cigar shape was 
basically represented by an octagonal cylinder. The rcunded 
nose of the torpedo was roughly approximated by reducing the 
diameter of the cylinder. The smaller diameter of a 
torpedo's tail was exaggerated by reducing the cylinder's 
diameter to a point. An approximation cf a propeller was 
attached to this point. As shown by the wire-frame image in 
Figure 33, even a crude approximation of such a complex 
surface required many polygons. The actual display cf this 
image required the definition of 27 polygons, 58 edges, and 
33 vertices. While the actual numerical values were not 
important, they shewed that the storage requirements 
increased rapidly with the complexity of the surfaces which 
were to be displayed. simple geometric figures, such as 
buildings, required little storage and the effort of 
approximation was minimal. The realistic approximation of 
complex surfaces required a large amount cf storage and 
quite alct of artistic talent. Finally, these numerical 
values indicated, as observed with nearly all 3-D figures, 
that the memory allocations required to store edge 
definitions were usually double that for any other image 
discriptor. 

The applications program, the polygonal penetration 
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subroutines, L00KE1 (the modified Looker subroutine) and 
their respective flow charts were included in Appendix B. 
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Figure 33 - TORPEDO APEROXIHI1ION 
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The Three-Dimensional Graphics Package contained in 
appendix A did not include two important aspects of computer 
generated graphics, an image shading procedure and 2-D and 
3-D graphs. Both items are important aspects of graphical 
presentations and should be included in a complete package. 
They were items left for future research due to the length 
cf time reguired to develop the present graphical software. 

The intent of this effort was to provide a portable 
computer graphics software package. To a very real extent, 
this goal was accomplished. However, the development 
software language, Fortran IV, was not entirely portable. 
In fact, a program which removed hidden lines on the IEM-360 
failed on the PDP-11 because of one its fortran 
idiosyncrasies. Additionally, an abnormal amount of time 
was expended attempting to input data via a file on each new 
computer. It would make more sense for an organization to 
really standardize the fortran supported by all of its main 
computers before a portable graphics software was 
implemented. 

Fortran as the development language for this graphics 
package proved to te quite efficient with cne exception. If 
Fortran IV had the binary operations common to languages 
like C and SAIL, the graphics software could have been 
simplified. While no time comparisons between two languages 
were attempted, the excellent, "real time", results acheived 
with the hidden surface algorithm indicated that Fortran was 
highly acceptable as the development language. Fortran may 
net be the best language available for graphics, but it is 


105 



the only universally supported and accepted software. 

The 3-D graphics package, as presented in Appendices A 
and B, provide a user with the ability to present any object 
on any selected display device which is supported through 
Fortran. Its portability was demonstrated cn four distinct 
types of display devices through the interface of three 
different host computers. There were only seven subroutines 
that contained statements which were device dependent. One 
of these, SHOffIN would not have been required if the raster 
scan CRT available had had square picture elements. These 
subroutines all call device procedures which performed the 
same task, but had different names. Bach of the device 
related subroutines caused the device to generate a line 
segment (cr vector) on the display surface. Usually, the 
display machines were built around a 16TBit display 
processor. The primary task of the device dependent 
procedures was a data conversion interface between the host 
and the display computer. Therefore, if a large 
organization utilized a "standard" Fortran and required that 
all interface subroutines utilize standard names, this 
Three-Dimensional Graphics Software could be made completely 
portable. The obvious advantages of this software 
installation would be a great reduction in software and 
personnel training costs. 
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APPENDIX A 


THREE-DIMENSIONAL GRAPHICS SUBROUTINES AND FLOW CHARTS 


Each subroutine was listed with a very brief explaration 
cf its function in the comment section proceeding the 
program. A specfic flow chart followed each program. When 
the contents of a program's flow chart could net be 
presented on a single page, a generalized flow chart was 
listed first. It was followed by specific charts which 
amplified all blocks of the flow chart that were marked with 
a circled number in a lower corner, such as: 

o 

The subroutines were divided into three groups. Display 
a Wire-Frame Image, Remove Hidden Lines, and Bidden 
Surfaces. The list of variables defined on pages 9 through 
16 were used throughout the programs and the flow charts. 

1■ Display a Wire- Fra me Image 

This group of subroutines included those to input 
the image data, the coordinate system transformations, the 
image clipping procedure, and the display subroutine for 
wire-frame images. Additionally, the interactive 
subroutines to rotate, scale, and translate an image were 
listed with their flow charts. The subroutines which were 
used to multiply matirces were included, but flow charts 
were not drawn due to their simplicity and lack of 
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CCCCCCeCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtCCCCCCCCCCCCGCCC 









c 



PARAMt' 
LINES, 
INI 1AL 


iMO.i r Ml MAGE All Hill 1 1 REMOVING HIDDEN LI NFS CALL 
INI UtSPLY A f rR M1I.U. utSPl.Y HAS NO CALllW 1 
MRS. 10 DISPLAY I HE IMAGE AElEtf PFMOVING ALL HIDDEN 
LALL SUOPOUlINt' KFMOvE (NO CALLING PAPAMLTLRS) AFTER 


1(1 USE 1H! SUtlROU 1 INFS ROTAU;, SCALI, AND THAMSL HIE FOLLU'/ING 
PARAMlRIO -NIST RE SPECIFIED VIA Tut USER'S PROGRAM: 

A. ROTAlt(Ii),IA*lS,Pl,P?,IHElA), ivHERf: 

1.1k- IS A IvM) ELE»UM ARRAY WHICH DEFINFS WHICH PANT 
01 I HE IMAGE IS ID HE ROTATED, AS FULLO.,3: 

f A 1 TO NO I »TE A SINGLE POLYHEDRON SE I ] K ( I ) - IN ( d ) - 
THE PULYHEOPON'S INDEX; 

CO) HI KUIATf SEVERAL CONSECUTIVELY INDEXED 
POLYHE DRAMS SF I : 

IN(I) = THE INDEX Dr (HE FIRST POLYHEDRON; 
irm) = the: index oe the last polyhedron; 

CO TO HOI ah the fmiipe image SET IR(I) = 31. 

?. I A* IS determines ruf a/is about .-.hum me polyhi 

OR ENTIRE 1 MAGE IS 10 ROTAIE, AS INDICATED: 

(A) 1 A X J S = 0 - ROTATE AHDUT OP-.JEC T x Axis 

(u) I AX ISa I - POT AIE ABOUT OBJECT T AXIS 

EC.) IA X1 S = P - ROTAIE ABOUT 0RJFC1 7 Axis 

ED) I A X I S - S - KOI Alt A fit,LIT AN ARBITRARY AXIS ■..Hir.i 

••'US! OF SPECIFIED USING THE ARRAYS PI AND IV. 

3. Pt(3) IS AN ARRAY WHICH CONTAINS ImE X, Y, AND 7 
EOiJJtCI) COORDINATE S (II UNI POINT ON ML ARiilRAHY AXIS 

EA) PI ( 1 ) IS ME x VALUE 

EB) PIER) IS ML Y VALUE. 

EC) P) (3) IS I HE. 2 VALUE. 

A. IVl3) IS AN ARRAY CONTAINING Ml X, Y, AND 2 
C(juki)i na ifs ue any nmm dieeeremt poimi un Tut. 

ARHTRARY AXIS, HUE PE ' V ( I ) , PREP), AND P<i(i) ARE USED 
FOR >., Y, AM) i f: Oi.'i-u I . A | E S, Rt'SPECrivELY. 

S. META IS ire: ANGLE IN DEGRLES THROUGH ivHICH Mir 
POLYHEDRON OP ME LMIRF IMAGE IS ID Hi. ROTAIFO, 

B. SCALE ( IN,S), '.MEWL: 

1. )R IS AN INTEGER ARRAY USED 111 DElcRMNf 'WHICH PAR [ 

UK ME IMAGE IS 10 HE SCALED AS DEE I NEE, ABOVE Mi- 


?. S (3) IS AM ARRAY v-. ITH: 

(A) SCI) IS Mil SCALE FACTOR FOR I He 

(B) SEP) IS THE' FCAl F FAC I OR FOR IHF 

EC) SEE) IS ME SEALE FAC TOR FOR ML 

C. I RADSL(JR,T), rtoFRF: 

I. IP IS Ail INTEGER ARRAY USED TO DEM.R■•'< 

Uf INF IMAGE io 10 TO. SCALED AS DF.EII.il) 

ROI A IE . 

?. 1(3) IS AN AnRAy '.-.I IH: 

(i>) Itl) IS EHF X DISTANCE TO TRANSLATE THF 
PUL YHE HR ON ' S 1 " A r. t . 

CEO T(<?) IS ME' Y Disr. TO TPf.'.SLATE ML 
POL rut 0RUN'S IMAGE. 

(C) 1(3) IS IhF 7 0 LSI • 10 TRANSLATE lot 
POL Y"t UPON ' S MAGI . 

THESE THREE SUHKOUT INFS MUSI OF CAlLfU AFIFR INI T A L AND HI F URL 
FIT UF H HISPI.Y OP REMOVE, INITIALLY. AF1IP 1 Hi IMAGE IS INPUT 
THE SJOROUIIMF INITAL NLFU .’it VIM Ilf CU.Irii AGAIN. MUS, IH 


■< coord mates! 

I. COURU! HAILS. 

! IF GtilCII PARI 
ABOVE Ml 









st- nuF mu 




w ,)i Tau i TTi N 



<E TE R ILOOK 

m'-'ICH IS USED in OE-ItR«JfIE WHICH SURFACES 
■ED AND WHICH ARE 10 HE REMOVED: 


I HI 


MIDI 


(FACES 

HIDDEN 


C WHICH ALCOwS 1 HE USER 1*0 SEE III: 

C 

CCCCCCCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

ccccccccccccrcccccccrcccccrcccccccocrcc a; ccccccccccccccccccccu: crtcccc 
c 

C 11-41 1 A L : INPUTS I HE IMAGE D A I A 

c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCriGCCCCCCCCCCLCCCCCCCCCCCCC.CCCCC 







I FOffGAI(qiq) 

? FORMA I ( JO I 0.5) 
J FORMA I (2I«> 
q FORMA I (Iq) 

S FORMA 1 tI0 10) 
h FORMA 1 (201 0 . J) 


7 FORMA!(11J 

e FORMA 1(?X, ' I HE KAUTEK DEVICE Vi OULU NOT OPEN') 
9 FORMA I (2X» 1 THE FUNCTION SCPkt'n FAILED') 

10 FORMA!(2X,'THE FU <t 1 ! 0*1 COl "•< T F AILED') 

END 


112 







































~ cccccccccccCtccccccccccc , c.cr.ccctcrcctcr.cCf.ccccccccccccecc(; CT: ccccccc<‘cc 


I Mi 


l C.L I PH 1NC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCC 
SUBROUTINE ROYCL.P 

INTEGE P PUIMTN,PClINTN*Ft*r,t h DG£M . t DOE ( 2 , 200 ) , f DGE 1 (100) 

INTEGER PCHANG (200) , FlJGt 2( t 00) 

COMMON /AAA/XE(l20),Yt(l?0),7MJ2U),P0INTN 
COMMON /AAH/ jfS ( I 2 0 ), V S ( 1 PO ) , IS ( 1 2 0 ) , Pf)XN TM 
COMMON /AC/ POOL 1,EDGE2, 1 ONE\ 

COMMON /A AO/ f UtsE /f.OGEM 
COMMON /CC/ Vji,VY,V/,A,LS,C/.CY,CZ 
COMMON /Ft /PCliA ,0, THf I A 
COMMON / J J / / C .< , V 3 Y , VCX, v C Y 

DIMENSION V (■', 0 ) , II (0,0 ) . I 2(0, 0) , 1 3 (A, A) , TO(0,0) ,RN((|,4) 

01 MENS 1 tl'i TMp.'K«),T£MPni) 

DATA V , T 1,T2,IT,! 

DO 30 1=1,120 


DO 100 1 = 1,LOOP M 
EDGECI,I)-tOOt1(1) 
EDGE (2,1) = i- 0GE2 (I ) 
160 CONTINUE 

00 13b 1=1,0 
T 1 ( I , I ) = 1 . 0 
13(1,1)=!.8 
10(1,11=1.0 
RN(I,1)=l.0 
136 CONTINOE 

T1 (0, 1)=-i/X 
T1(0,2)=-VY 
T1(0,3)=-V7 
T 2 ( 1 , 1 ) =- 1.0 
T2(2,3)=-l.0 
T2(3,?)=l.0 
12 ( 0 , 0)=!.0 
SO=YX*VXtVY>VY 



T0(2,2)=CSS 
To (2,3) = -SSi; 
TO(3,2)=SSS 
To(3,i)=Cl)S 


sx=i .o 

SY= I .0 

IF(VSX.GT.VSY) a<=/sY/VSX 
IP(vsx.lt.vsy) gy = v:.x/vsy 


■ IN- 

















































cc£ccccccccccccccccccccccLcca;ccf: , rrcC£rccc.nccu cccccccccccccrccccL^ctt 


C CLIP: CLIPS T HP IMAGE AGAINST THE VIS 1 RU tljSPLAY 

C AND ELIMINATES ALL PORTIONS OF EDGES <HICH are off 1 HE 

C SCREEN 

C 

CCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCLCCCCCCCCGLCCCCCCCCCCCLCCCCCLCLCl 



COMMOh /AAB/ XS( 120 >,YSCI 20 ),ZS(I 20 ),P 01 NT': 

INTEGtP POLYGNton, II I , ^ DOE ( d , 2 0 :j ) , PC M AMG ( 2 ' 0 ) , PI'LGN > EDCEM 
S,SHAD( 60 ) , PUR. TP 
DU ISH 3 = 1 / EDGE 



J( 1 ) =f.DGE ( 1 , I ) 

J(2)=FDGfctf,l) 

J1=J(1 ) 

J2 = .) (?) 

X ( ) T = >: S C J 1 ; 

X(2)=XS(J?) 

Y ( 1) = Yt> ( J 1 ) 

Y (2 ) = Y S ( J ? ) 

Z(1)=7b(J1) 

Z(2)=7SCJ2) 

I TOT(1JrJCUDL(X(1),t),Z(1),I ) 

HOT(2 ) = I CODE(X(2)/Y(?)/Z(2),2) 

148 I 0= T T 0 I (1)11 TOT (2) 

IFdrt.lO.OJ GO 10 1'U) 

DO 111 K = !,4 

lArJCMMl ,K)tICHK(?,K) 

IEdA.FU.2) GO 10 14? 

1 A I CONTINUE 

JF tlTOTCiS) . f.AI. 0) CALL 1 Si- AP (IS, L N) 

JR=J(IS) 

PCHAflG (1 R) = 1 

IF (ICHM IS, 1 ) .tO.O) GO TO III A 

TS(7(IS)EX(IS))/((X(1SJ-X(LM))-(2(LN)-Z(lS)) ) 
Z(1S) = I*(?(LN)-?(I5)MZ(1S) 

X (I S ) = - 7 ( IS ) 

Y(1S)=I*(Y(LN)-Y(1S)) + V(IS) 

GO TO 14? 

144 IF (TCHK OS, ?) .KJ.O ) i,0 TO 1 US 

T=(7()S)~X(IS))/((X CL'U-MIS))-(ZCLN)-Z(IS))) 
Z(IS) = 1*(Z(I IW(JS)H/(1SJ 
X(IS)=Z(IS) 

YfISJrf *(Y(LN)-Y(IS))♦Y t1S) 

GO TO 14/ 

1 us if(ilhi\us, 3 ).cd.O) (,n to iu« 

I = (7(IS)<Y(IS)}X(O(IS)-Y(l N))-(Z(LN)-Z(IS))) 
Z ( !S) = [ * (7 (I. N)-7 ( IS) )< / (IS) 

Y(IS)--Z(1S) 

X(l3J?l>(X(tlS)-*m,))U(IS) 

GO TO 1U 7 

1 4 6 T = (7 ( IS)-Y ( IS) ;/ ( (Y (LII)-Y (IS) T-(ZfLN)-/ ( IS) ) ) 
ZtIS)*»*<Z(LN)-ZtIS))NIIS) 

Y(1S)=Z(IS) 

x (ir.) = t * (x(un-x( is)) m (is) 

14 7 iiui(i5) = itnuC(x<is), y n; 


), 7 (IS), i: 
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ccccccccccccccccccccccccccccccccccccccccrcr.cccccccccccc(.ccccccccccrLro 

c 

C ICODE; USED fir CLIP TO DETERMINE. IF rtN tND P01M IS VISIBLE AND 

C IF NOT TO r.HICH SIDE OF I HE DISPLAY THE POINT IS LOCATLO. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCL 
FUNCTION 1CIJDI TX.Y.Z, i 1 
COMMON /II/1CHK(?,9) 

00 159 J=1 ,n 

159 1CHK(I,J)=0 — 

1FCX.LT.-Z) ICHK(I,n = l 
TFCX.GT.Z) I CHK t T . 2) = I 
IF(Y.LT.-Z) 1CH*CI,5) = I 
IF(Y.GI.Z) l CHK ( I , '4 ) = I 

I CODE-T CHK f 1, 1 ) * TCHM \ , d) ♦ I CM". (I , 3) 11 CHK (1,9) 

RETURN 

END 
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ccccccccccrcuf.u-LC( rr.r.rc Lrtrccu:crf;c:crccc :crccrccccccccccccccccccrt 

C IS *ttP: S'-.APS IHl VALUES Of TWO JtiTFJfHS. 

C 

CCCCCCCCCCCCCflCCCCCCCCCiCCCCCCCCCCCtCCCCCCCwCCTCCCCCCCClCCCCCCCCCCCCC 

I T = I 

J = I T 
RtlURM 

two 


ccccccctccrcccctci.rcctrcrcctctrccLcccc.c.f:rf cccccc t.LCCcccccLccccccf.Lcc 

c swap: soaps thf values of iwo real numbers. 

c 

ccccccccrcccr.t a cccccccccteeccccf Grccctccccccccf.ccccEccGCcccccccct.cG 
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c- 


IT=J 

J =1 

I =IT 

^ RE r 

MS A! 

T=X 

X=Y 

Y=X 

rURN ^ 

(X,Y)^ 

^RET.UI 

3N ^ 
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ccccccccccccctcu.cccccrccccccr.ccccccrcrtcccccccrcccr.cccccccccccccLCccc 

c 

C SCRM : TRANSFORMS EYF. COORD I jjA T FS TO SCREEN COORD I N A IE5 

C 

ccccccccccttcccc.ccccccccctccccccccccccccccrcccccccceccccccccccccccru.c 

SUBROUTINE SCRN 

COMMON / A AH/ XSC 1 ZO) , YSC 1 20) , ZS( 1 <?()) .PfllNTM 
COMMON /FF/ Pf.HAMO 
COMMON /JJ/VSX,VSY,VCX,VCY 
INTEGER PCHANUC20U),PCH.\HM 

IEtS.GI.VSY) S = VSY 
00 IS I 1 = 1 .POtUl ' 

IF(PCHANGt J).CO.1) GO 10 ISt 
Z=ZS(I) 

XSCI)='S(1)/Z*S«YCX 
YS(I)=YS(1)//*S«VCY 

zsi n=-i .o/z 

151 CON I 1 i'll It 
RETURN 




















ccccccccccccccccccctctccccccccccccrtccccccccccccccccccccccccccccccccr.c 


C 0I5PL.Y: MASTER SIlURUUIINE WHICH DRAWS THE CLIPPED IMAGE (IN I HE 

C SELECTED UUTPUI DEVICF. 

C 

CCCCCCCCCCCCCCCCCCCCCr.C'CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SURROUT INE UISPLYUK) 

DIMENSION IRf?) 

COMMON /AA/POLYliE , POLYHN 
COMMON /AB/f-OL Ylilj, ”l')lbN, SHAD 
COMMON /0 AD/ EDGC , ERGE M 

COMMON /AAH/ XS(IPO)i(StliO),ES(120)» PO1 NT M 
COMMON /FE / EDGCMl, (200) 

INTEGER POL YPN ( «>0, 1 1 ) , E DGE (2, 20 0) , POLPN , EOGEM, COl JR 
INTEGER VEC I OH , SHAD (60 ), PO I N ! M, F.ObCHP f PUI. YI1E ( 1 0 ,»? > , POLYHN 
CALL RDYPLP 
CALL CLIP 
CALL SCRN 
DO JO I - I,EDGEM 
JO EUGCHGin-0 
JL-IR(I) 

if( jL.ra.in go to is? 

KL = 1R(?) 

T SPOLYHE(JL,t ) 

J = POL YHtIKL,?) 

GO 10 IS?? 

IS? 1=1 

J-POLG.N 

15?? JLMzSHAD(l) 

11 =COl.OR (JLM) 

DO 1SJ0 K=T,J 
LN = POLYGN(K,I 1) 

DO I5A0 N=I,LN 

L = POl YGN IK, NT 
KlrEDGE(1,L) 

K2=f DGE ( ? i L ) 

lEIM.EU.Ol PO TO IS'lU 

IF (EDGCHGd ) .EO. () PO TO 15-40 

EDGC HG(L) = 1 

XSI=XS(M) 

YS1=YS(M)/?.0 
XS?=XS(K?) 

YS?=YS(K2)/?.0 

KP=VECTOP(>St,Y5I,TS2,YS?) 

IF(MT.LI.O) ivR I I E l 6, J ) 

1SA0 CON fIMUE 
I 5JO CONIINUF 
RETURN 

J FORMA 1('I HE FUNCTION VECTOR FAILEU'I 
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ccctcccccccccccccccccccccccccccccrcccccccccccLrccccccccccccccccccccccc 


01 ATE : R01A1ES t 


COMMON /At!/ POLYGON 
COM,•ION /AC/ EuGtt.Ei 
COMMON /AAA/Xt(120), 

COMMON /TF/PCHACC, IHI-1A 
DIMENSION 1 ( U , A ) , I-! 1 l A , « ), R21 <1, A ) 
&LMC«), LR C« ) , P I ( 3) , P2! 3) , 1 R(2) 
INTEGtR POLY HE C 1 0, 2 ) , POL YC.'. 10 0 , ! 
(tl*OL rHri,POLGM,lO:.E N.PGINlfi, PU1 11 ,1 
DAI A I,k|,R2,R3,11, ILMP/Rb‘0.0/ 
DO 30 1 1 = 1 f P01\ I Is* 

PCHANG(111=0 

1HElA = lof *3.H1M/1«D. 

JL=1R(1) 

IF(JL.E0.31)GO 10 111 


GO 10 [112,113,HA).RAX1S 
do ns k = i r, j 


call KOlXltoGHCL)) 
CALL ROIXCEUGE2CL)) 
llh CONTINUE 
115 CONIINUE 
RETURN 

112 DO 118 K=11,J 


CALL R01Y(FDGE1(L)) 
CALL kOT Y(EDGE? CL)) 
IIP CONTINUE 
110 CONI I fill!. 

Rill URN 

113 DU 120 K = I T , J 

LN = PUL.YCN(k, 1 1 1 
DO 121 N-I.Lh 

L Id = POL Y GM ( X , M J 
CALL R()1?(( L'GL1 (l ) 1 
CALL kill 7(EDGk?(L ) ) 



11 A DX=P2(1J-PI CI ) 

DY = E‘2l?)«Pl (2) 

D/ = P?C3)-Pt ( -,) 
DIFMP:DX*OXtOY*OYin/«L)2 

F.X-DX/D1SI 
EY = UY/I)1ST 
E?sO?/niSl 
V = SGRI (I Y*FY<L7«( L) 

DO 122 L = 1 . 

1 (L,i) = I .n 
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H !l L , L 1 = r:B 
R?(l , L ) = 1.0 
R3(L,L)= 



T(«#1)=tM t I) 

1(4,P)=PIIP) 

T (iJ, 3 ) = P1 (5) 

R 1 (P,2)=«A 
R1 (2,3)=-l<p 
R1 (3,<M=»il 
R1Ci,i)=HA 
R?(l, n-Rh 
R2(1,3)--R0 
RP(3, I )=KG 

R2(3,3)=RP ^ ^ ^ 

CALL GMPUut I P.i-.P, T , II , ■/, A, A) 

It (I. (0.31) GO TO 123 
DO I PA K =11,J 

00 IPS 0-1,1 fi 

L=POLVO ifK.M) 

CALL k(> t AWHU OGt 1 (L ), T1 ) 
CAIL ROlARiitrOOtPtL) , IT ) 
IPS coon out 

IP4 CON) I L’Uf. 

RPI UR m 

! ? 3 DO 1P(. H=1,I'NIM(N 
CALL HOTAWIUk, Tl ) 

1P6 CON II fUH. 

111 oo io n,' h, ipR, I pi) , 

DO 1?7 k=»,P01NTM 
CALL KCilMM 
j ? i cun non 

Rf TUwr. 


rf'AXlS 









128 DO 130 K-l.POHITM 

CALL ROTf(K) 

130 CON1INUE 
Rfc TURN 

129 DO 131 K=1 ,(*niNTN 

CALL RuTl(K) 

131 CONTINUE 
RtTURN 
END 
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131 




































32 




















































33 



















































cccccc'c'cccccccccccccccccccccccccccrcccccrr.ccxcctccccccccccc'cfcrcct.ct n 
c 

C ROLL: USED HY ROTX, RUTY, RuTZ, AND ROTARR TO PERFORM A RfcPfcAlfl 

C CALCULATION, THE COSINE AM} SINE OF THE ANGLE OF ROT AT IOI.. 

C 

ccccccccccccccccccccrccccccccccccccccccccccccccccccccccccccccccccccm 

FUNCTION ROLL(T,P,AN) 

ROLL=I *COS (AN) t S « 5 IN (A) 

RETURN 

ENO 




ti0LL=T*C05 (AN )+3* JIN(AN) 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccl 

c 

C ROTX: PERFORMS THE REPEATED PUl TI PLICATIONS T 0 ROTATE THE 

C SELECTED IMAGE AHOuT THE X AXIS 

C 

CCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCnCL 
SlIHROUTJNL ROTXCLR) 

COMMON /AAA/ XE ( I d\) ) , YE ( 1 ? 0) , ZE( 1 20) , P01NIN 
COMMON /EF / F'CHANG, THETA 
INI EGER PC HANOT <? 0 0 ) , POINT A; 

IECPCHANGTLRT.nl.) T PL TORN 
Y A-YE C L R) 

ZA = Zt CLP) 

ZE(LRJSROLLC/ArYA,I hL IA) 

ZA = -Z A 

YLCLR)=R'QLl (YA,ZA, TmETA) 

PC HANG Tl N) = l 

RETURN 

END 
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ccctrLCccccLcccccccLn cc:ccLCccicc.a:a:( cr.Lr.crcrccLcccca c rcLCLrccc-ccr. 
c 

C ROTY: PERFORMS THfc RE PEA TER RULI.'S 10 ROT A 1E THE SELECTED 

C IMAGE AHOUl THE Y AZI3 

cccccccccccccccccccccrccccccccctccccccrcccccccciccccccccccccccccccccc 

SORRCM.IT 1 ME JffiTYlLRT 

COMMON /AAA/ xE(l<?0),YE(l<?n),ZEU20),PUlNTN 

COMMON /FT/ PCHA-.G, THETA 

INTEGFP PCHANG .POINT'D 

IF ((’CHANG (IR) .Ell. 1 ) RETURN 

XAsXEU.fi) 

ZA=ZE(LR) 

ZE(L R)= ROl l (ZA,xA, THETA) 

ZAs-ZA 

XE(LR)=ROLL(XA,ZA,THETA) 

PCHANGTLPJsI 

RETURN 

END 
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cccccccLccrcccccccctcEcCTTncCccccccccefdfccLCCctfCccccccccccrcccctrcrc. 

c 

C ROTZ: PFREORMS IHt REPEATm MULl.'S in ROTATE THE SELECTFU 

C IMAGE ABOUT IHt ^ AXIS 
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iccccccGcccccccccccccecccccceccccccccccccceccctcccccccccr.ccccccti 


SCALE: scales / 


,LE POL y rifcCR:)-‘. ( 


CCCLCCC.CCCCCCLCiiXCCCCCCCCCCCCr.LCLCCC.LCC 
. 3 C*Lttia,s) 
lUfi S( 3 ), Xkt-’) 



:cccctr,cccccr.ccccLCCCCCccr.ccc 


, PCHA.'Jt (<?U0) 


te»ULYGfi(*«A»» 

CALL SCL (f-'Mc 1 tL ) J 
CALL SCLdOCF^fD) 


133 COM I Ufc 
RETUPIj 

I it DO l 3b A=1 , i->ri 1 >j'I fj 
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ccccccccccccccccccccrcccccccccccrcccccccccccccccctcccc.ccccccccccr.ict.cc 


C TRANSL: TRANSLATES A SINGLE PCLYHFDRON Ok I HE ENURE IMAGE 
C 

ccccccccccccccccccccccccccccccccccccccccccrcrccccccccccccccccccccccccc 

SUBROUTINE l'RANSL (IR , I ) 

DIMENSION T ( 3 ) , I R ( ? ) 

INTEGER PULYHt t 1 0 , ? ) , PtlL YGN ( 6 0 , 1 1 ) , EDGE 1 ( 1 0 0 ) < F DGLP t 100), 
fcPCHANG (<)0U) , POL'.’., EDGE N, Ik AS I 1 , POIN1 N, PULYHN,SHAD C60) 

COMMON /AA/POt YMt» POLY‘• 't 
COMMON /AI) /P 01. Y GN,P 0 L G N,S H A 0 
COMMON /AC/ f-l’Gt! , EDGE 2 »! ’'GET-. 

COMMON /AAA/ xf ( lEO), Yl. ( I <?0) ,/E ( IPO) .POINEN 
COMMON /FE/ MChAAG,IPLI A 
COMMON /EL/) %, 1Y,1 7 
7 X = T ( 1 ) 

1 Y = T(2) 


00 30 11 = 1, P0 ININ 
30 PCHAMGCII)=0 
JL=IR(I) 

IEtJL.ro.3)) GO TO 103 
KL=IK(?) 

J=P01.YHEtKL,P) 

I I=POLYHK(JL, I) 

DO 105 K=l,J 

LN = POLYGNtA,I 1 ) 

00 106 M =I,LN 

L=POLYGNfh,N) 

CALL NE x 7 P T tF DOt I tL) ) 
CALL NEXlPIttOCLPU. )) 

106 CONTINUE 
105 CONIINUE 

kE TURN 

103 DO 107 A = 1,POINTN 
CALL NEXTPTfK) 

107 COM 1 I .Milt 
RETURN 
END 
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TRANS (IR ,T jl 






ccccccccccccccccccccccccccccccrcccccccccccccccccccccctcccccccccccccccc 


C NEXIPI: PERFORMS REPEATED OPERATIONS FOP TRANSl 

C 

ccccccccctccccctccctrcccccccccccEcccrcccccrccccccccccccccxcccccccctcct. 

SUBROUTINE uEXTPTILR) 

INTEGER PEHANC-(POO) jPOINTN 

COMMON /AAA/XE(120),TE(1 PC),ZE(120),POINTN 

COMMON /FF/PCi'ANG , 1 hF. I A 

COMMON /Ft/IX,Tr,IZ 

]F(PCHANG(1KJ.FU.1) RETURN 

XE(LR)=XE(LR)-TX 

YE(LR) = rF.UP)-TT . 

ZEII.R)=ZE(LR)-TZ 

PCMANG(LW)=1 

RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccrccccccccccccccccccLcccL 


ccctcccccccrccCLCcci 


crcccccccccccrLccccci.rcccccccccrcrtccccrccccLCi: 


cccccccccccccLCLCi:c(.rccLCCccrcci.cu:a;ccccccccr.cccccurc;ccccccccrc:cLn 


I 10 1=1 

I) = o. o 
I 10 J=l 
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2. Hidden Line Re moval 


Included in this group were the subroutines utlized 
to remove hidden lines from the wire-frame figure produced 
by DISPLY. The subroutine REHOVE initiated this procedure 
when called. 




ccccccccccccccccccccccccctctccctcccccecccccccccf-ccccrccccccccccccccrct 

c 

C REMOVE: MASTER SUOROUT!W AlUCH DKArtS THE IMAGE, AFTER REMOVING 

C ALL IIIOOEN' LIMES, ON THE SLLEC TFO LfllTPUT nF.Vil'E. 

C 

cccGCcrcccctccccrcccccrcccccrcr.ccLccccrcL'cccrcr.cr.ccccccnccr.cccrcrLrccc 



IRDYCLP 
SCRN J 
lim'kkm(m) 


LI 

Hi* 


EPS=0.001 
IXX=VSX* 2 + EPS 
IYY=VSY* 2+EPS 
1=0 
J=0 


V ARNCK(I,J,IXX,IY' 
DISPL2 


■55 
















ccccccceccccccccEcccctccccccccctccecccccc^ccccrcfccecccccccccccccccccc 








EDLINK(K)=0 
TERMINATES LIST 
SORT EDLINKO 
LIST SO THAT: 


POLLNK(I)=POLPTR 
POLPTR=I 
J=P0LYGN(I,11) 
LR=POLYGN(I,1) 


IFNDEDCLP 

~r~ 





^ RETURN 
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:6o 




OLDP=J 
GO TO 215 
213 PULPTR=J 
GO 10 210 
212 OLOP=P 

P=PULLNK(P) 

GO TO 215 

211 If-(CHANGE.EU. I) GO 10 203 
RtTURN 
END 
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INITI.-.(POLGN) 




ORDER POLYGONS IN 
POLLNKO WITH CLOSEST 
POLYGONS TO VIEWER 
LISTED FIRST 


OLDP=P I 
P=POLLNK(P) 
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sill 


cccctra 


: rarctccTracrarcccececracracmTrrarrrcrc c trcrrcc cm 


c 

c 

c 

c 


WARMCK: lb [Ht 
ARE HIDDEN. 


MAIN 


INE WHICH BE TERM] 


[Nf.S «H1 


fCH LIMES 


ccccccccrcccccr.Lrtcccctr.ccccccccrcrcrcccccccrcrcccccccccrcccrcccti.CLCL 



LS17EX,ISI7F 


common /rij/ puli '•*> polp ir 


COMMON /RC / POLY A C t»u ) , PUL Y H CbO ) , POL YC ( 6 0 ) , POL YO (60 1 , POL 7 MM ( of ) 
COMMON /RO/ XP, yP, ZP, xm, YH, Z\ 

COMMON /RE/ POLL SI 



339 WLX-Ltn 

SI?F=ISIZFX 

WRX = feLX*Sl/E-tPull.M 

SIZEMSIZtY 

WuY^lriOI T 

WTY = r'HY*SIZt-fPSILM 


t N T E R = 0 



333 Ir(P.tQ.n) GO TO 318 

IF (POl ZMNIP) .(ir.ZMli'.MX) GO 10 318 
the:ta=o 

NEXTEL'-POLFBGCP) 

331 CALL GFtNEXU) 

IF( I . t 0 . 0 ) GO TO 819 
CALL CUP 8 CJR) 

IF(JR.FO.O) GO 10 330 
POLLSI(P)=INTLH 
1 NUR = P 
1 HETA =-1 
GO HI 319 

330 IHP'lA-THFTAtOEl IA 1 
GO 10 331 
319 CONTI Nllt. 

IF ( 1 AilS (THE f A ) .E 0.8 ) GO 1 U 838 
P=POLLNh(P) 

GO 10 383 
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ZMAxl: 


ZMAx2 = 0• 
7P.AX 3 = 0. 
7MAXU=0. 


227 ir (SUKRND .E 0.01 RtJ 10 ?2'J 

Z2 = GEIZ(S‘JR>no,*LX,.'.TTl 
Z3 = GEIZFSUP..\0,s.KX,aSF> 

Z« = GE I Z (SIMP'Jl),P x , i, 1 Y ) 

if- r/1 .g r ,/y r. ii go to 225 

IFC/Z.GE.Zf-'i-V) go fo 22S 

lF(Z3.r.t.zv[M) go rn ?25 
JF (Zil.Gt.Z’AiPNl c.0 10 22S 


Z MI M 1 = 7 1 
ZMAX1 -7 1 
Zf-;iN2 = Z2 
ZMAX2=72 
Z M l hi 3 r 7 5 
Z M A X 3 = 7 3 
7M1N9=Z9 
?MAXU -Zu 


GO 10 226 

225 if (zi .LF.z«#xn nu in 297 

IF (Z?.LE.Z u Ax21 GO 10 792 
IFF Z 3. LF..Z x3) GO H) 2«2 
1M/'1*U.Z“S*'0 GI 10 292 

226 SURPMU = pnLL5nSURi<\D) 

GO 10 227 

292 PENE T=I 

IF (Z 1 .L 1 f.'l JZ'-lNl -, 1 . I 
IF (7 1 .F.l ./MAxnz"A’Cts/ I 
IKZ2.U .ZMi*'iY')/l- , l*!2 = /7 
IFC/?.GT.ZX'*X2 )Z-aX 2 = 7? 
IFIZ3.I. I . 7 P I f‘ 3 1 /. '■ 1M3-Z3 
If F/3»«'1.Z«A> 3 )Z.«mX 3=/3 
IF (Z<f.LT .Z^r-'A ) Z v IN'l=i<l 
IF(Z 9.G1.ZPAXu)Z v Axa = Zu 
GO 10 226 

229 If fPFNEI .f,0.1 ) GO 10 226 


I IFfP.EO.OJ GO TO 229 
II fll]G‘K.r.U.07 i.O TO t 
Z I =GE I Z (P »L X , •. r V ) 

Z 2 = GF 1 7 (P« XL *»I V ) 

7i=Gl 1 7 (•’, A><» ,a*i ) 

Z9 = GRIZtP»ftHx,r.lY) 
IKZl.Lt./MAxn CO in 
IFIZ2.LL ./*’AX2) GO ro 
IF (/ 5.1 F . /Y'U U l.n 10 
if fZ6.lt ./MAX'- I GO 10 
J = P0I.LS1 CP) 

IFF ill UP .t*). 0 1 GO TO 231 


.r.F .ZM.INI1 
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1?2 







































m 


ICX=(ISIZEX) MODULO 2 
ICY=(ISIZEY) MODULO 2 
ISIZEX=ISIZEX/2 


ISIZEX=ISIZEX+ICX 


//5HOWIT 

y(WLX,WBY, ;lx,wby, 
4shad,d / 


ISTPTR=I5TPTR+1 
JJ=ISTPTR 
ISTACK(JJ,1)=LEFT 
ISTACK(JJ,2)=IBOTT 
IST/ CK(JJ,3)sXSXZSX 
ISTaCK(JJ,4)=ISIZEY 


WHEN 1=1: 

IBOTT=IBOTT+I5IZEY 
ISIZ SY a IS1Z Y-IC Y 


ISIZEX=ISIZEX—ICX 


I=ISTPTR 

LEFT=ISTACK(I,1) 

I BOTT=ISTACK(1,2) 
ISIZEX=ISIZr;x( I ,3) 
ISIZEYalSXZEY(1,4) 
ISTPTR=ISTPTR+1 


'■/HEN 1 = 3: 
ISIZEY=I.iIZEY—ICY 
ibott =1 [a tt-i.;iz;-;y 















































CCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCtCCCCCCCCCCClCLCL 


c getz: determines the depth of thf selected f-olygon at any puint 

C ON thf display screen, i.e. thl sreen coordinate z value. 

c 

CCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCC 
FUNCTION GEIZ(I,X,Y) 

COMMON /RC/ POLYA (60), PUL YOK'D), POL YC (60), POL YD (60), POL ZMN(oO) 
GETZ - (-POLYA ( 1 ) <-X-POLYP ( T ) » Y-POL YO ( I ) ) /POL YC (I) 

return 

END 
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(getzCi,x,y)^ 




GETZ=-[pOLYA(I)*X+POLYB(I)»Y + 
POLYD(I)J / POLYC(I) 




^ RETURN ^ 
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ccCFccccrccTcccccccccCCCCCCCCCCCT'CCCCCCCCCftcccccccccccccccccccccLnrt. 

c 

C CLIP2: CLIPS tHt SF.LFCTED EUGt AGAINST 1 ME CURREN1 DISPLAY 

C WINDOW IF IMF. FOSE INTERSECTS THE WINDOW OR COMPUTES TmF 

C ANGLE SUhTENDFD BY THE EOGfc WHICH IS USED TO DETERMINE IF THE 

C CURRENT POLYGON SURROUNDS THIS WINDOW. 

C 

CCCCCCCCCCCCCCCCrCCCCCrCCCCLCCCCCCCCCCCCCLCr.CCCCCCCCCCCCCCCCCCCCCCCLCL 
SUHHOUI1NE CLlPfCl) 

COMMON /RG/ DELTA I , wl. X , y.R X , SPY , .IT Y 
CUMMON /RD/XP.YP, /P, XB, YN, 7N 
COMMON /RTI/ XX(2),YY(2) 

COMMON /PI/ ICHK(2fR) 

COMMON /P J / IC(2),!S,1N 
INTEGER DEL TAI 
DELTA I=U 
I iM=2 
1 S= I 

XX(1)=XN 
XX(2)= XP 
Y Y ( I ) = Y N 
YY(2) = YP 

icc i )-jcodf(xx( n ,yy cn,n 

IC(2)=JC00t<Xx(2), YYC2),2) 

2S7 icoi = ir(i)•ir(2) 

IF t ICOT.EU.0) GU TO 2S0 
DO 55 KK=1 , <4 

I AAA = )CHK ( I ,KK ) ♦ ICHw (2»I\K ) 

IF (1A A A. FI 0.2) GO TO 251 
55 CON 1 I NOE 

IF (IC(IS).FN.O) CALL I SNAP(IS,IN) 

IFCICIWT IS, 1 ) .ET4.U) GO TO 253 
CALL PUSH(0,WL X ) 

GU TO 25 7 

253 IFHCHK(1S,21.Eu.u) GO TO 250 
CALL PUSHTO, V,RX) 

GU TO 257 

250 IF(ICHK(IS,3).Ett.O) GO TO 255 
CALL PUSttT 1 ,WisY) 

GO TO 257 

255 CALL PUSH U, WT Y) 

GU TO 257 
250 I = -l 
RETURN 

25! I A 1 = I ANGLE (XN, YI41 

1A2=I ANGLECxX(IS),YY(IS)) 

I A 3 = 1 ANGl L (XT’, YP) 

IA1=IA?-]A! 

IA2=I A 3-!A2 
CALL IDEUIAI) 

CALL I 0 E T(I A 2) 

DELTA I:I A 141A2 
1 = 0 

RE TURN 
END 
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CCCtCCQCCCCUCCCCCCCCGCCCCCCeCCCCCCCCCCCfCrCCCCCr.CCCCCCCCCCCCCCCr.irLC 

C IUEf: USEO riY CLIH? TO CDHHICT I HE SUHTEN0EO ANGLE OF AN tDl.F. 

CCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCCClC 
SUBROUTINE IDEItl) 

IF(IABStn.LE.S) RETURN 
IFd.GE.O) GO 1.0 2oO 



260 I=1-p 

RE TURN 
END 
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CCCCCCCCCCCLCCCCCCCCCCCLCLClr. 


' rcccccccccccrcccrccccrccccccccccccccccci 

C JCODE: 10 Otltf-VI.VE IE 40 r\r PHIM of AN FPr.fc IS VISIBLE 
C IN The CURhF .1 .Moon.-, AO,, if NOT l'("-'PUIES THE POSITION 10 IHL 

C END POINTS (K. THfc SCPFE'!. 

C 

EECEECCCCCCCCCECCCCrrLCLCCrCrCECCCECrCCCrCCCEOEEECCCeECECLEEECEOCLTLEC 

COUPON /PL/ l.n LI A I , „LX , . PH , ¥. !i T , ATT 




































cccccccccccccccccccLCcrccccicccrccrcccccccrccccccccLCCccccccccccccctcr. 


ccccccccccccccccciccrcccccccrccccLccccccrcccccccccccct .ccccccccc can ct 

SUBROOTI ML PU8H(1A,BT) 

COMMON /Pi./ X<C^> r VM (21 
COMMON /P.I/ ICC2) ,IS, IN 
IF ( I A.FJJ. I) Go IP 27S 

Yy(iS) = (YY(iM)-yyns))«cHi-xx(is))/(xx(in)-xxus))tYYiis) 

XX( I S ) = H T 


275 ) 


27b 


S) = (> 


(C IM)-XX(IS) ) » CHI - YY US) )/(YY( IM-YYllSn tXX ( IS1 


yy i is)=m 

276 IC(IS) = JC.U1)t(XXtIS),YYnS).lS) 
RE TURN 
FND 
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ccccccccccccccccccccccccrcccrcfccccrccccccrcccrccccccccLCCccccccccc'cc 

c 

C SHOW 1 1 : STOWES lHt IMAGE in the APPROIAlf ARRAYS FOR DISPLAY 

C BY SUBROUTINE UISPLYJ. 

C 

eecccccccccccccccLCCCCccccccctccccccccccccrcrcrccccccccccceeccceeccLC 

SUBROUT I fife SHllwJTlXI , Y 1 , it, Y j, I , J ) 

COMMON /DA/ XSS(?000),YSS(?000) 

COMMON /DB/ JO 

JO=JDt1 

XSS(JO)=X1 

YSS t JD)=Yl 

JD=JDt1 

XSS(JO )=X J 

YSSUOJrYJ 

RETURN 

END 






^SHOy,'IT(Xl,Yl t X2,Y2,I,J) ^ 


JD=JD+1 
XSS(JD)=X1 
YSS(JD)=Y1 


JD=JD+1 
XSS(JD)=X2 
YSS(JD)=Y2 


^ RETURN ^ 
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■cccccccctrcci 


cceccccccccLccccccccccccccccccccrcccccrcccrcrcrtcccccLf 

C DI SPLe > : ACfuALLY ORA „s THE I»AGF NITh The hUHTIN U 

C UN THE SELECTED DEVICE. 

C 

c c cc ccctccctccrcccccrtTCctf.LCLCC r c c cc c cc cc c c c c ( :cc E cccccrcc 

CUMMfJt'. /DA/ XdSl?00V) , Yi>3(?00v) 

COMMON /DH/ Jt; 


DtS RtMuvtO 

ccccccccrccc 


INTEGER VICTOR,COLOR 
DATA tPS. ACC'/0.I'OO 1 ,«?00.0/ 
J=CUI OWl ISHaD) 

DO AOU 1=1,JD,? 

11 = H I 


X1=X5S(T ) 

xa=xssn i) 

Y!=YSSM )//.(> 


T I ON i/EC TOR T Mil 
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RETURN 





3. Hidden Surfaces 


By calling the subroutine SURFACE, this group was 
utilized to display the object with its hidden surfaces 
removed, or tc display the hidden or back surfaces, as 
determined by the calling parameter ILOOK. 
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pru =FU Gur.t 




5b IF(PTR.fQ.O) GO TO 55 
NFXl=l UILS1(PIR) 
lrP(],P|R) 

IF(J.NE.O.ASb.SHAUd).NE.O) GO 10 57 
I=P(?,PIP) 

IF (l.i«F..O.AwD.SF'ADU).F‘E.O) go TO R7 


57 


IF ( YS< J ) .Lf .YSl 


I-Y3FJ) tO.< 
1 F ( I . L I . I . ( 
ENT LSI f *' T rf: 
I Yfc'vTK f ] ) =f 
PTR=MfcX( 

GO TO 5-, 


i T .TYPES) GO TO 100 

i t v n j 


IF CSFG.f-O.O) |,0 TO (.0 

XLE.FT (5tG) = XLLF I (oF>.) HxIlF I fGFG) 

XRIuHT f„; G) H)Xk.,l! \ ( SI G 1 

7Ltr 1 ( SLG J =71 ■' F r fSFJi.) i 071 t M FSf G) 

?kluHI (afi,) = 7PlGil fSF-vi) tC/PGHl fjFG) 

lYl-lYLtf T toFG)« I 

1Y2-IYPGHI(GFG)tl 

I YLL F I (,SCt>) = T Y 1 

I TRGHI C Fife 1,1- I I? 

TK IY I .Gt .0.4 :0. T rP.Ofc.U) GO 10 01 
PlR-POLGIJ.JClGi.) 

IF (PTiY.LO.OJ GO TO o?. 

IFCCHAF.oh(RIR).GK.O) GO TO 61 
CMAlvGt (i 1 T k ) -POL CHU 
PULCHGsPIH 


CALL RFX.SKT (St.G) 
CALL HE I H L X (Si U I 
SEG=IxSKGT(M G) 
GU 10 b« 


RIP:IYt KTH(iY 1 
II (PTK.LCi.O) GO If) t,b 

Ivvi =f dgf;c i ,pir) 

I VV r>= t nilF (i>,P!R) 

ITT lib I =YR fWVI) 

Iylas i = rs(i nvi 

II.'ELY-I Y F -< S I - I rLA.'jf 
RL ALOY = Y S ( 1 i/Vc')-YS( I V V 1 ) 
IF f I Otl.Y .I.F.O) GO 10 „? 


XSL UPL = (>S(itfV?)-AS tIVtflIJ/HE ALOY 
XF 1 RSI = <5( I v V 1 HXSLOPL MY-YSCIVVll) 
2 SLOPE-CZGtlVVP)-iS(1VYI))YRtALOY 
7F I RSI - 2S ( 1 V V 11 + 7 GLOPt * l Y -Y f, ( 1 VV 1 ) ) 



IF (PI . F. u . 01 GO 10 6M 
CHAOIH (PI UPOLCKG 


StG = .SLGLS I (P I ) 






IF(SEG.LQ.O) GO TO 70 
IT E 1-0 

ir ( (XT IRST .Lt .XLf FT (St GI) .OH. ( (XF I RSI .EO.XLEtT(SEG) I 
(XSLOPF .1.1 .OXLEF T (StG) )) ) I It I = E 1 
1TF?=U 

IF ( (XF I PS I . L F. XP I GFM ( SEC ) ) .OH . { ( XF 1 ‘ S I . 10 . XR I GHT ( St M 
(XSLOPF .LI .OXkGFiT(SFG)))) I T E = £l 
I Y I =0 

TFC1YLFFT(SEG).LI.0)Ir1=? 

I Y?-0 

IF (1 YHGi'T l SKGJ .LT.0 jIYS = 1 
ii =j miYM 

GO in 5,; , 3, 70, S, 'U. 5,3, 70,7P, 5, 70, 70,70) , I I 


SEG=PUt St G t SI l 1 
GO 10 71 
SLGl- IGIHLM I • I 
If CSEG1 .EO.O) GO 10 Io0 
POlGO'UStGI )=P,I 
X Lt F1CSfcGl) = »F: PSI 
DXLt.F I (SEG I ) = xi LOPE 
7LEFT(SEG1)-7F 1 k S I 
DZLE F I (SEG 1 1=7. ILOPE 
I YLLF f (St G1 )-! Iff Y 
CALL PUlXST(St U1 
IF (PPEVIS.F.0.01 GO TO 75 
PlTLStG(PktVJS) = tf GI 
GO 10 7« 

SLGLSi (PI ) = StGI 
POLSF G(Sf 0 1 ) =SE 5 
GO 10 <>» 

SEGl = 1GIbLK(TP 1 

IF (StGl.EO.O) GO TO 160 

POLGO F(SFGI ) = PI 

XLt F !(St G 1 1= x LF F ’(SEG) 

OXLEF T ( SEG 1 I = OXL o I (SLG) 

7LF.I T (S( GI )=7LEFTi SLG 1 

07LLF f (SFGI )=07LtF 1 (St I-1 

I Yl.LF- I (SFGI J = 1 Y L£ ■ T (GLG) 

IYI EF I (SFOl^O 

XW1 GFIT (SFGI JSXKIR *V 

DxRGMI fSFGl 1 = XSLOPl 

7FMGHT (Sl l.l )=/r IPSI 

tlZKGFII (of GI >=/SLOP- 

I YIHitFI (Sf GI ) = 1 OF I.Y 

CALL l'OIXSF(StGl) 

ir (PPt VIS.1.0.0) Go in mill 

PULSt u(ppr.viS)-.sF>.i 

GO TO I.U'I 

SLGISKPI J = S F G I 


CUM 11. ■ J l 
PTP=ENTLSI(kIP) 
GU fU b', 

PRLVIS=StG 


51 G = P 0LSt I, (St (.) 


1 GO TO / r j 








PRF vi: 


StG=StGLSl(PI) 

7b IFCSEG.tCJ.O) GO IU 66 

Irl-IYtfcFTCSEG) 

IY2=IYRGhICSEG) 

IF(IYI.LT.0.4ND.IY2.LI.0) GO 10 95 
IF(IY».EO.O.AIiD.IY?.EO.O) GO .10 77 

iffm.we.e.OK.jY?.GE.o) go n m 

IYI.EF7 (5E G) = IY PUH F ( 3FG) 

n«GHT (SEG)-O 

XLE.FT CSFlOlixkifOfl (St K ) 

DxLtF I (StG)=l)U>GHI (SFG) 

7LFH (SLG J = /P I Gh I (StG) 

D71EF1 f St G} -D7 r GH F (S! G ) 

GO JO 7o 

77 1 SPOUSEG«StG) 

IF (PPttfJS.Fu.y) GO 10 79 

POLS[G(HkcVISJ=i 

GO 10 PO 

79 Sl.GLS I (017 = 1 

80 CALL RP/IxT (Sf.G ) 

CAt L PLKsLK(SLG) 

GO 10 7u 

78 Nf X l--POLStG(SfcG) 

IF fi'.FxT .E u.yj CO 10 11(1 
IFCIYLtFT (l> t x l ) .Gl.01 i;-J 10 hi 
IYHGH 1 (t»E (> 7 = 1 YLfcf 1 (NEXT) 

IYLtFlCNF x T) = 0 

XK1GM1 CoF o)=XLFF T (AEx I ) 

DXKGFI1 f SF G ) = ITxLt F I C.-.l XT ) 

7RIGHIf&FG)=/LEFI(NfcXI) 

D7RGHI (SE. U) -O/LfcF 1CNEX1 ) 

GO m 76 

81 IF flYt-GnllM x I ) .GE.O) GO 1(1 82 
lYfiGHI (SI U ) = I Y..GH1 (.'IF.xT ) 

XKJGHT fSf.i,) = XWH,HI (:. F X T ) 

DXflGOT ( 88 U ) = OXPGH I f fjF XT ) 

7WIGH1C S t G) = 7 k1G H1(NEXT) 

D2PGH1 (oFG)=lV.RG'iT (NEXT) 

GO ID 76 

82 POlSEG(SFG)=P)LSEGCNf xT ) 

CALL W"XSkT (N't* I ) 

CALL RfcfKLuCHtXI) 

GO 10 7o 
7b CHG = 0 


86 IF (LEG. L 0.0) GO 10 85 

I = 1 xsfg i r sEi.) 

IF (I .80.0) 00 10 8 5 

IFfxLLF1 (oCG).Lt.XLFF 1 ( 1 ) ) GO TU HS 


K=IXSLF1(StG) 

IF(K.M.'.O) 1 Y 8 a 0 I ( K ) = 1 
1 XSLF I ( n-K 

I XSLF | (if ( r )=l 

IF (K.nF.0)1«SIFI ( KI - St G 
IxSHGI (bti.)-K 
1XSKGIC1)=StG 
IF ( SF l, c o I . F N . SEG I 


SFGFRIrl 








0 7S 



GO III <V 

iMxirr irsno.r.r.xsf'NRG) on to iom 
CURSF (.= I 1 lil'G l iCUHSrt.) 


1.96 



f (AltF I (t,tG) .L 1 . I .0) GU If t 
F ( X L K F I (f>l G) .GI .XHLSL) t.U U 

FfiPU.NE.LSTSFG) GU TO 600 


CALL i'F'XSR I ( SFG ) 
CALL RF. I i;LK ( SLR) 


x a x = xi_eF i (si r,)-* I .o 



IF f X IF X (in 1 rst 0).GF.XSS) GO TO (>01 
IAC I VL(St'») = SU-uoT 
IF (SFGOUl .F.J.0) SEGLO-Sf.G 
SEGOUI-GfG 

GO 10 (Rvo,«07I,1L0OK 
CALL LoOixF.K 1 St G ) 

GO TO GO 

CALL L ."’Mil fSHI) 

GO 10 Oil 

I ALT VE ( 3t.fi) = SF.*»A‘GT 


CALL UllNiLPl 1 IR,SEG) 
GO 10 (i.Oi, 160) , I 10 
IMbkGOUT.F.I.O) GO TO 
IACIVLC>FGLil)=SLGACT 



Fl.tQ.O) 

Sill C I»’PLF I ) 

ALL Kt I FU.K 11MPLF T I 
IFIXSPNHG.LT.XRFSL) I 
T.E'I.O) GO TO 


TO < 


IF ( 


GO 10 611 
610 SAMF5T = (V 

fall CALL SHOW 

S9 CONIlM.lt 

Rt T OWN 

100 rill 1 T F 16# 1 d I ) 

1?1 F0RM1 (’F.KPJR: IMF Y COORUINAIE Hi U F F I ME SCREEN.') 
no «h 11 r.» 6,111 j 

1 1 1 FORMA I 1 ‘ERROR: Nl XT = 0 ■ ) 


I 60 Y.R1 IF.U.W1 
99‘> FoRMAI Id*. 1 I Ml 
ILL INI • ,/fX, • lilt 


NI.IMHLR OF SrORtMS GENfRAT FD FOR 
STORAGE PROVIDED.’) 


IMIS SCAM 
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r SURF AC \ 
UlSR.ILOOKM 

CONVERT OBJECT COORDINATE: 
TO SCREEN COORDINATE > ^ 

ZERO ARRAY IYFNTRO, LINK 
' ) 

AMLNK CN THEMSELVES 


9. 


ORDER ACTIVE LIST OF 
SEGMENTS BY INCREASING 
X VrtLUE OF THE LEFT 
ENDPOINT 


2 ) 


SEGCNT=3F.GACT=0 
AND REMOVE ALL EXITING 
IMPLIED EDGES 


2 ) 


I INITIALIZE POINTERS 

EXAMINE EDGE AND 
DETERMINE SCAN LINE EACH 
ENTERS THE DISPLAY. STORE 
THIS INFO. IN IYENTRO AN 
ENTLSTO ( 4 


IY=I -IY IS THE SCAN LINE 

POINTER- 


UPDATE SEGMENT ENDPOINT: 
USING SLOPES OF THE X & Z 
•D LINE EQN.'S AND MARK 
IE POLYGON OF AN EXITIN< 
SEGMENT AS CHANGING FOt 
ACTIVE SEGMENTS 


COMPUTE SLOPES OF THE X ft 
Z 2-D LINES AND THE INTER¬ 
SECTION OF ALL ENT RING 
EDGES MTH SCAN LINE IY 


PLAC ; EACH ENTERING EDGE 

01 CHE ' C : /'a 

SEGMENT LISTS _1 7 / 


INITIALIZE POINTS AND X 
VALUE OF FIRST SPAN 
(XSPNRG) TO ZERO. USING 
PRUVIS SCAN LINES SAMPLE 
POINTS DIVIDE SCAN LINEJ^ 
INTO SPANS. /% / 


PUT ALL SEGMENTS WHICH 
INTERSECT THIS SPAN INTO] 
TV'O LISTS, SEGOUT AND 
SPGArT- ^ 

HAVE LOOKER OR L00KR1 

EXAMINE ALL SEGMENTS IN 
THESE TWO LISTS TO 
DETERMINE BOX COUNT SrtTDV. 
RCX TYPE _(j-3 J 

HAVE THINKER DETERMINE 
MLin' IF THE CONTENTS OF THIS 
SPAN CAN BE DISPLAYED,—. 



REMOVE ALL EXITING EDG :S 
FROM POLYGONAL SEGMENT 
LISTS, AND CONSOLIDATE 
SEGMENT Bi. S' - r ! -GS ( 


LOCK SI 

75 
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PTR=0 


fNO 
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|V . 
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ccccccCGCcccccctr.i.ccccn.cccLCi:ccccrcccc( ctccccccctcccccccccccccccuC 


C lGTeLn: ]b IJStU TO Olll/Ufi THE T'inLX'Of A FREE SI KMENT OLOLK 

C OF STORAGE. 

C 

CCCCCCCCCCCCCf CUCCCCtCCCCCLCCCGCCCCCCCCCCCCCCf CCC.CICCCC.CCCCCCCCCCCI. 
FUNG [ 1 ITT'i IGIHLKU) 

CO'' I'-O.M /I A / I AC r Vfc (ijO) , IF RE LG 
ClW'On /In/ TYLLF I (SO), lYRGHT (bill 
1-1T RtLd 
IGTbLA=U 

If(i.fcO.n) HE 1 URN 
1 YI.EF I (I )=l> 

IYRGH! f J T = 0 
TFRtl S= l M; 1 vF ( I I 
1G T U L. A = 1 
RE I URN 
END 
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RETURN 






Q RSTBLK(I)) 




IACTIVE(I)=IFRELS 

IFRELS=I 




RETURN 


216 









If |-§ 


cccccccccccccccLccccccccrcccccccccccccccn.cccLr:cccccccrcccr.cccccccrT 


c 

C COLOR OR SHADE ON THF RAMTEK SCREEN. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCL 



INI EGER StGCM , POL SfeG, JHILGON, CPI VOL , PUL UN , SHAO , vEC T OR 
INTEGER COIOP 
ISAI-tPsO 
y i = i y 

IFISFGCHT .L I . IT RETURN 
DO 701) IM.SEOCNI 
Xl-TSAf R 
ISEG = 1SSEGI I ) 

X2=Jx 

IFUSEO.E. . 0) GO TO MO 
IP = PGLGniiUSEG) 

JLK = SHADUP) 

KL^VECTORUl, yi»x?,yit 
iF(KL.LT.O) GO TO ISS'l 
710 ISAMRrlX 
700 CONTINUE 
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318 






















ccccccccccr.ccccLcicccircccccccctr.cLCcu.cccr cccccccccccccccr.ccccu i 
)), rsc ivo j ,^s( 1 eo) , nnnmi 
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rCtttccccLCCccr.cci.CLa.rc ccctcccCCCCCCcctcc ccccrcccccccccccccccccctcc 
c 

C LOOKED: EXAMINES Int ACllvL L 1 S I OF SEG-lE HIS AND DETE DP 1 ML S 

C HUft MANY APE V 11 .< AKLE IN TmI 3 SPAM. IT- ALSO OF NEK AIES 

C » HO* AROUND ALL VIE.-iAHLE SEGMENTS AMD PASSES THE 

c NUMbC«< of segments ape in the hox and hox type iu the 

c think*. 

c 

cceccccccccLCccccccLa.rcr.cccccccccrLCCccrcccccrccccccLCccc'cccccLf;LCL 

SUKWOLMINt LOOKER!M 




ZLFf I ( fe 0 ) .ZRJGHI CoO) 


m/KGHI , HZ'" I <1 . it/'- Ax 
S/PNmI 

V, loFHl L, ISFULL 


IF (S.UFF T.GI.xSPNLF) 
S/Li.Fl=Z|r,T UGPNLF 
SXLtFI = *SPMl F 

no ID 6 I i 

611 ISFUl L = U 

S*Pt.MT=Z Ml <X6 !! ': ! 'G) 

sxrghi =xsp'n<n 
no io si.'i 
613 I SFIjLl = 0 
6 1 A SDTV = SXLF.FI 

IF CSXLEF T.l L.XSPMLF) spiv-: 
if (iy xcfjf.ml.d 5 no rn sis 

CALL l OmIIiFXI I ) 


GO TO !>1 1 


GO 10 ols 


6 1 S 


If (JHXCMT.Mt . 1 ) CO in 616 
Z 1 = H Z i N I (SXLFF 1 ) 

If ( (IUIJ F I ,lF .S*Lt'F I ) . A,o>. filxPOMl ,r, t . .SYRGHT ) . AND. 
8 (Z 1 .LE.SZI.FF T ) .AND. (Zc'.lF.SZPC.FII ) ) KMUHN 

Z i=ZIHlUHAl l-F I ) 

?? = Z IMF (FURGhI ) 

If CSXLEF I 1 .•'XLF'F 1 J no TO 616 
if (SXknnt .L 1.0VKG1I1 ) no ro oi« 

IF fZl .C.I . 'i7L F f 1 ) nil ID ol(> 

IF (Z?.ni .K/hOMI) 1,11 Til o 1 6 
CALL LliAhuXCIJ 


IF (131 DLL .I'D.D) Nil 10 o|0 
J F C I HI ULL .> i:.») Go HI 6H> 

T E w P = liXL E f T t (IJXRGiil-i'.L EM) MS7LE FT-HZI.f FI )/ (dZDGHI - 
H/.t EFT -SZKGitl 1 3ZL F. F I ) 

CALL XPNiWtt) 

im=lt:"P 

II (liZLFt r A F.S/Li.F 1) CALL I itr.M' ( 1 tlSLG I , I USE G? ) 

K1 TIJPM 

CALL *PiJOll*(l) 


IF (SXLII T .01 .11*11 F FJ 1.0 10 6?0 
IF (SXHGItl .LT .hXMF.Ill ) i,0 TO 6?il 


?21 


616 
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ccccccccccccccccci:cr.ciccci.ccrt;ccru.tcccc.ri:ccccccccr.cr.ccccLCc;ccrcci((. 



SZU F IS/1.LF 1(1) 
S^KGHW-Mr.nr t t) 


ISKULL=J 

IF (bXLfcF 1 .G I . XSPt.LF) (.() 10 (j 11 
S^LLF I =n -IK XbPULF ) 

SaL LFI=X5PHLF 
GU 10 hi? 

Ml ISFuLL=0 

b\? n (SXPGrfl .L I .XSP.JPG) GO 1U u 1 5 
S/PGH I =/ 1 ill UbPNPii) 
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-SS 





IF(SZLfcf T .LI.o/MAX) GO in 620 
IF (SZHGMr .LI .II7MAX) (ill 10 hd 0 
CALL LOA0i>X ( I ) 

RhTu»w 

IF (13VLEF I.GI.SxLfcFT) CD lU o?l 
IF (bXRGHT • L I •oXkGMT) GO lu 021 
IF (tiZMIN.L 7 .SzLtFl) GJ 111 621 
IF (tiZ^TN.LT .SZKGHT ) GO 10 621 
RFTURM 

CALL xPNI)bX(l) 

RE TURN 
END 
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ccccccccccccccccccceccccrccccccctccccr.rcrtcctcfcccr.cccccccccrca.c4.tc 


C PUTXST: PLACES THt SEGMENT RASPED TCI II AS \HE CALLING 

C PARAF'EILR I INIO UF THE XSOR1 USI3, IXSLFT(I) 

c and ixstiurcn. 

c 

ccccccccccctCLCCccccctrcccccccctrccccwoccccceccccccccccccccctirtccec 

SUBROUTINE PUIXMU) 

COMMON / TC / 1 XSl.F T tht)} , I X3RG f ( 60 ) , SFGFST 
INTEGER SfcGFSl 

iF'EStGFST .ME . U) 1XSLF I (SE SFST ) = I 
TX5LFI (l) = f) 

TXSKGT ( 1 ) =SLC.F SI 

SC5FSI-I 

RETURN 

END 
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ccccccccccccccrccccrrccccccrccccrcccccccrcccrcccccccccccccccccci 


C IMF VI'KMULt btQMLNIS I U tHt CllKFFtfll SPAN. 

c 

ccccccctcccccccccccr.cccLcccLCCcLCcccccccerer.ctcccca;ccf.cccccr.crcCLC 



/TF/IFSXCNT# Ii)Xl YF> 

/tC/SXUM .8XRGH I 
/TH/SXLFFr,S*iiGMT,S7LF.FT ,?ZBUH[ 

7 T 1 7 I tf:Sf IjJ , raster,oil/,3U1V, P>F jll, ISFIJLL 


IBXCf«T = II>3iCHTt 1 
IKSD1V.LI . u Iv),)lV = 
IF (SXLliFT .L ! ,'iXIJ F I 
IFCGXF'GMT.C-I ..iXRC-T 
IFFSZLF.F T .L I .H7-M [.'.'1 
IF .LI ,f-7>'Vl) 

IF ( S7LF.F T ,G r.d/Vfx) 
IF(S/i'Ghr.GI.iJ7rtAX ) 
RLTUOH 


END 


) h X L F F .t =S XLF.F I 
J SfXRGH I =Sx»GHT 
IU -IjiM-S/Lf F T 

F57..:i.xr<5^1 FF I 
H7FA x = 37KGF* F 
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cccclcccclclclccclclccccccclcccccl ccc'ccl-Clccc c l ccccci.ccc l rcc l r l rt i 

c 

C THfNKR: RECEIVES IMF. INFORMATION ABOUT fHL CONTENTS OF 

C THIS SPAN hHILH NAS 0£ VU.OPI.P BY THfc LOOKER ANO 

C DETtRMlNES IF THIS SPAN CONTAINS 0] SPLAYAliLF. 

C INFORMATION. 


ccccccccccrcccccr.ccccccccr.cccccccccccccccLccccccccccccrccLccccr.cr.LC 

COMMON /I ■ J /I Yl.FFT CSV) ,IYRGHT CoO) 

COMMIT) /T C/1X SLF I (On) , IxSRiil CoO) , SF.OFSI 

COMMON /TD/XLfcFT (60) ,XRK;mT (60) ,7LEF T (60) ,2RIGH1 CbO) 

COMMON / 11 / x SPNL E , X:iPN>'G, xRFSL 


I /TF/IiJaCGT, 11' X I YE 
i /TC./ilXL EFT ,4 xPGH 1 ,l(?Lf 
i /I I X 1 f<j[ ij 1 , luSLGcJ.OlV, 
i xtj/i str.sM. r--.Pi i , impi 
l / Ik /1 PL St. G( 60 ) t tPLGONI 
XIL/'JXlT E I (b'l) ,0*RGH I ( 


P/RGHT , i)7 -IN,ei/nax 
V, I> E 01 L, ISHJLL 
I P«l vS 


IF ( lBXCNI.Nt .0) CO 10 0?3 
CALL STUPlPCXSPNWG,1?LP0) 


6?3 if nnxCNJ. i) co m bin 

IF (0 xlF.F 1 . Nt . X SP'TL F ) L A LI. S 1 OKF P C t* X L F f T » Ut KO J 
CALL STUPtPCfUKGi.I , 10SF.G1 ) 

IF (OXNGHT .Nt .xSP'iNG) LALL G I ORF P (XSPW:, J r<0 > 
CALL RFCOkOC T*'Sf.C 1 , iONF , loCfc I ) 


6?6 


bib 

640 


6 31 
6?B 


1 F (I fix T yP. f.v . I ) co 10 6<?s 
CALL S T 0 ■< t. P ( P t V , IhStGl I 
CALL WF.CfWI) (I USE G ) > IDNEj (ZERO) 
ISE0S,M=i'-SI C 1 * I UOOO * 1 iiSr CP 

1 - IMPLSI 
I PPL VS-0 

IF(I.F0.0) GO TO bin 
IF ( I St. GSM. FU. I PL St G( D) GO 10 bib 
IPRf VS=1 
T=IXSNGIfI) 

GO 10 6? 7 

IF ( I .1 0.0) GO 10 nib 

if (iPwrvSito.cn no in m'» 

IXSWGI (1P«EVS)s(Xo»i»I C I ) 

GO 10 6.30 
I MPLS I = 1 XSNit T i ] ) 
oxi.EFi m=r> iv-xleei c u 
xLrri c n=i>iv 

OXX-Xl.feF I ( I HuXLFF T ( 1 ) 

1 r (OXx.L 1. 1 ) c,u I il o 41 
IMPXx .0 I . /Pi SL ) 1.0 III 0 3 1 
CALL P'lIXSI C I) 

CALL l-FCOWO ( I, JOIN, J 71 Wo) 

GO 10 6 4? 

CALL PFIRLK(I) 

GO 10 63? 

J »IGTHl'K (III) 

IF(J.LO.O) 00 TO 160 
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CCCCCCccccccecctrccccLCtccccccctccccccccccccccccrcCLCCCCCtccCLCLCtrc 


ccccccccccccrra.rrnccLr.tc 




636 


II-( 1 .nr. 11 ) i:.j 
X = Xlt'F I C ISI’LF 
CALL PUfSA-’U 
I MPLF T = U 
11 (lYLtFl[I ).61 


1)00 TO t.3? 



call pursA'n*) 
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ecccccccccccccccccccccccccccctcccccccccccccccccccccccccGcceccccccccc 


C STORFP: USFU BY THF THII/KR TO RECORD 1 HE DISPLAY DATA 

C FOR THE CIIPkTkT SL'A.'i L I MIt . THE XS VALUES ARE 

C STORED IN THE ARRAY TSPUSUJ AMO THE SEGMENT'S 

C INDEX IS STORES IN 1S3TGC1). 

C 

ccccccccccccccccccccccrcciccccccctccccccccccrccrcEccccccecccccccccrc 

SUBROUTINE STOREPTX,!) 

COMMON /TO/SEGCM.L STSEG 
COMMON /TP/ISPDSTi.O) , ISSEGCtU)) 

INTEGER SEGCNT 

IFC (SEGCNT .NF.O) .AMU. I I-EU.LSTSEG) ) GO 1(1 610 
SFGCNT = SF.GCNT < 1 
LSTSFG=I 

640 ISPOS(StGCNIJ-X 

ISSEGTSEGCNI1=LSTSEG 

RE TURN 

END 


Zb 9 
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APPENDIX B 


APPLICATIONS PROGRAM AND SUBROUTINES 


Ihe applications program and its flow chart were 
presented first. Next, the seven subroutines which enabled 
the user tc determine polygonal penetration, change the 
coordinate values of a polygon, and alter the shade cr color 
of a polygon were listed. These programs and their flow 
charts were not included witt the 3-D graphics software 
package since their applicability was strictly related to a 
tracking presentation. 
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ccccccrcccr.ccLCCcccLccccccccctccct rcrccccLrccccccccccccccccccccccLCL 
c 

C THIS IS AN APPLICAIiON PROGRAM DESIGNED TO DEMONSTRATE THE 

C USAGE OT THE THREE-DIMENSIONAL CONPUTtR GRAPHICS SOFTWARE 
C PACKAGE. THE PROGRAM SIMULATES THE TRACKING OF A TORPEDO IN AN 
C IRREGULAR TORPFOO TESI AREA. 

C 

CCCCCCCCCCCCCCCLCCCCCLCCCCCCCCCCrCl Cf CfLCCCCtT.rCCCCLCLrrCLCCCCCCCLCC 
DIMENSION I K ( P) , S ( 3) , I (US), T t 5 1 , T w ALK X f S ) , THACKY (S ) . TR ACK L ( b ! , 
&P?f 3) ,PVECH)( 5],MT),PHU 
DATA UVFCTCi/P. 0,-1 .0,0.0/ 

CALL INITAt 
CONV=lftu.»/3. I'URR 


1 = 1 

IH {c*) =5 


INITIALIZING I HI. OLD MKHIJOU VECTOR 


IT IS VISIBLE 


IRfI)=# 

IR(^)='4 

P2(1)- 10 . U 

P?(P)=10.0 

P?(3)=I0.0 

CALL SCALE(IR,PP) 

ROTATE THE VfHlCLt SO THAI |T ]s INI TALLY PROADSIDE TO THE 
VIEaF.R. 


THElA=R0.0 
I AVIS-? 

CALL ROI A1F (IR, I Ax IS, PI ,P<?, (HI TA) 
REA0(S,1) T 
DO 10 .1=1,3 
10 StJ)=-T(J) 

1 AX IS-3 


BEGIN THE RECURSIVE PORIIG* OF I HE PROGRAM 


POO TRACKX(|)=!tlI 
tracky m = i c?) 

1 RACK Z(I 1 = 1(3) 

DO PO J=1,3 
P(> T Q (J ) = T (J) 

C D1SPI AY T HE lURPEP'i 1 tS 1 ARE A KITH HIDDEN SURFACES SHOWN. 

C 

I LOOK -? 

IR (<?) = 3 

CALL SURF AC l IR. It. OOk ) 


1ANSI.AIE 
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ccccccccccccccccrcccrcccrccccccccccccccccccccccccccui rcccccccrccu t 


C CHGSGN: CHARGES THE SIGN OF THK COEFFICIENTS FOR Hit 

C PULYGUN K. 

C 

ccccccccccccccccccccrircccccrcccccccccccccccccrcccccccccccccccn rt (i. 

SUBROUT I NF CMC.SU'J(X) 

CO^MOi'. /CA/ PUL YA (f,0 ) , POLYP (60) , PCIl.YC COD) , PULYD(hO) 

POL Y A( K)--P01 f A ( K ) 

PlILYBIK )--POLYTUK j 
POL YC(K)r-PULYC(Aj 
POL Y 0 t K J = - PlJl. Yf) CK ) 

RETURN 

END 
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cccccccccGcceccccccECCf cctctccccccrccrtccc.rcrcccci'cicf.ccccccci'ct C( 


C INTERN: COPRIMES IHE POINT OHICH IS UiR GE ORE T PIC CENTER 

C OF A Sfl OF rurivrx POIYUEURA. TrlEf), IT CHANGES TOE 

C SIGN OF THE COEFFICIENTS FOR ALL PTU YGOrtS OF EACH 

C POLYHEDRON SO THAI I HI OO T PRODUCT OF AN IN TEH I UR 

C PUINI ivJTH THi VECTOR NORMAL WILL Bl LESS THAN ZERO. 

C 

cccccccccccccccccccccccccccccccccf ccccccrrccor.cc r.cccccccccccccctrc a. 


01 Ml 


, /CA/ POL YA(6 0 5 
/C.u/ *SI , YSU'- 


I, Pol YCC60), 


EDGE ? 160) ,FDI 


< E DGK2,EOGLN , 


1220 

1210 


KL=Ik(2) 

DO 1200 I=JL,fcl 
ML = E'OL YuF C 1,1) 

Nl rPOLYHF(1,2) 

YSUN=0.0 
I C M = 0 

lo=polyC'*: ia , i I) 
DO 1220 N-|,L N 


IE m-nANGCK I ) .E.,.0 J CALL SU-AMIKKI) 
IF (PCHAND iK2 1 .1 0.0) CALL gUM'tTlK?) 


CONI (MIL 
C=ICN1 
XAVE = xSUR-/C 
YAVE=YSOH/C 







INTERR(IR) 














CCCECCCCCCCtCCC«UCrcCCCtCCCCCt.E**XftCCtXfcCCCCtCCCC:CtCCCtCl.CeCCCtCt 



COMMON /AC/ EDGE I (I Ufl) ,LDGt/( 1001 ,E06tN 
COMMON /AAA/ X! 1120), YE 11/0),/f;n?nj .P0IN1N 
INTEGER POL you , Priix-'J, f 0GE1 ,tOGEP, toot N.POIM IN 
IKJ.tO. I I )oO Hi 211 
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jit:um(ip, j,p; 


■=POLYGN(IP,J) 
I=EDGE1(L) 
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1 Phong, Bui Tuong, "Illumination for Cuirputer 

Generated Pictures." Communications of the ACH, v. 16. no. 
6, p.311, June 1975. - 

2 Schmacker, Robert A., Sproull, Robert F., and 

Sutherlsnd, Ivan E.. "A Cnaraterization of Ten 

Hidden-Surface Alqorithms," Computinq Surveys, v. 6. no. 1. 
p.9, March 1974. -*-- 

3 Kewman, William A. and Sprcull, Robert F., 
Principles of Interactive Computer G ra ph ic s, McGraw-Hill, p. 
4 5B—252, "197i. 

♦ Ibid., p. 247. 

s Ibid., p. 256-259. 

* Ibid., p.249. 

f Ibid., p.316-317. 
a Ibid., p. 545. 

« Scbumacker, p. 13-14. 
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